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ABSTRACT 

Based on the ultraviolet to far-infrared photometry already compiled and presented in a companion 
paper (Barro et al. 2011a, Paper I), we present a detailed SED analysis of nearly 80,000 IRAC 
3.6+4.5 ^m-selected galaxies in the Extended Groth Strip. We estimate photometric redshifts, stellar 
masses, and star formation rates separately for each galaxy in this large sample. The catalog includes 
76,936 sources with [3.6]<23.75 (85% completeness level of the IRAC survey) over 0.48 deg 2 . The 
typical photometric redshift accuracy is Az/(1 + z)=0.034, with a catastrophic outlier fraction of 
just 2%. We quantify the systematics introduced by the use of different stellar population synthesis 
libraries and IMFs in the calculation of stellar masses. We find systematic offsets ranging from 0.1 
to 0.4 dex, with a typical scatter of 0.3 dex. We also provide UV- and IR-based SFRs for all sample 
galaxies, based on several sets of dust emission templates and SFR indicators. We evaluate the 
systematic differences and goodness of the different SFR estimations using the deep FIDEL 70 /mi 
data available in the EGS. Typical random uncertainties of the IR-bases SFRs are a factor of two, 
with non- negligible systematic effects at z>1.5 observed when only MIPS 24 jxm data is available. All 
data products (SEDs, postage stamps from imaging data, and different estimations of the photometric 
redshifts, stellar masses, and SFRs of each galaxy) described in this and the companion paper are 
publicly available, and they can be accessed through our the web-interface utility Rainbow-navigator. 
Subject headings: galaxies: starburst — galaxies: photometry — galaxies: high-redshift — infrared: 
galaxies. 
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1. INTRODUCTION 

Multi-band catalogs are the fuel for studies aimed 
at exploring the global evolution of galaxies over cos- 
mic history. They have been used to study the red- 
shift evol ution of the star formation rate (SFR) den- 
sity (e.g.. iHopkins fe Beacoml 120061 iReddv et all I2008L 
iBouwens et al.1 I2009D. and the stellar mass ass embly 
process (e.g.. Bundv et al.l 12006 1 iFontana est al.1 12006L 
iPerez-Gonzalez et al.ll2008l iMarchesini et aT1l2009h . 

The unprecedented sensitivity of modern surveys de- 
tect millions of distant galaxies to faint flux levels that 
for all practical purposes lie well beyond the capabili- 
ties of even the most recent multi-object spectrographs 
at the largest telescopes. As a consequence, their in- 
trinsic properties must be estimated through multi-band 
photometric data using fitting techniques to stellar pop- 
ulation templates, and/or empirical relations. Among 
the basic parameters needed to characterize a galaxy, ar- 
guably the most important is the redshift, which must 
be inferred from an analysis of its Spectral Energy Dis- 
tribution (SED). Photometric redshift techniques are 
now sufficiently accurate to derive statistically reliable 
conclusions for high-redshift galaxy populations (e.g., 
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iSilva et al.lfl998l iWolf et al.ll200l lllbert et al.l[2009h . 

Many different codes have been developed to calculate 
photometric redshifts based in the same principle: find- 
ing the galaxy spectral template best fitting the observed 
photometry in se veral band-passes. Som e exa mples in- 
clude HYPERZ (jBolzonella et al.ll2000f ). B PZ (iBem'tezj 
2000h or LePHARE (Arnouts&Ilbert; e.g., lllbert et al.1 
2009). The implementation is very sensitive to the qual- 



ity of the photometry and the capability of the observed 
bands to probe key continuum features of the spectra 
(e.g., the Lyman and Balmcr breaks). It also depends 
strongly on the availability of templates that are sta- 
tistically representative and successful in characterizing 
the emission of galaxies. The impact of these factors 
in the uncertainty of the estimations is not straightfor- 
ward, and it can lead to catastrophic erro rs beyond the 
simp l e propagation of the stati stical errors (jOvaizu et al.1 
2008, Hil debrandt et al.ll2008D . In recent years, several 
techniques have been developed to improve the relia- 
bility of_the photo metric redshifts (e.g., Bayesian pri- 
ors (iBemted l2000f) . template-optimization procedures 
(lllbert et al.l l2006a[) and mach ine-learning neural net- 
works (|Collister fc Lahavll2004l )V Recent work including 
some of these advances have achi eved remarkable p reci- 
sion [e.g., Az/(1 + z) < 0.012 in lllbert et al.l 120091 and 
Az/(1 + z) = 0.06 at z>1.5 in Ivan Dokkum et alj|20"09j . 

Once a galaxy's redshift has been estimated, the most 
significant physical properties that can then be derived 
from multi-wavelength photometry are the stellar mass 
and the SFR. However, estimates derived from mod- 
eling of the observed SEDs involve significant random 
and systematic uncertainties. The estimate of the stellar 
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mass by fitting stellar pop ulation synthesis models is a 
wides pread technique (e.g . jBell et aL 2003, P anter et al.l 
20071 IWalchereTa l. 2008) that requires making some as- 
sumptions regarding the initial stellar mass function, the 
star formation history or the extinction law. Moreover, 
there exist significant differences among stellar popula- 
tion libraries. These differences can lead to discrepan- 
cies in the stellar ma s s estimation o f a factor of a few 
(|Maraston et al.ll2006l [Bruzualll2007t) . 

SFR estimates based on UV and/or IR luminosities are 
considered reasonably robust for large galaxy samples 
with multi-wavelength photometry, wher e other trac- 
ers, such as spectroscopy, are unavailable ( Reddv et al.l 
2006, iSalim et al.l [20071 IDaddi et al.l 120071) . A major 
problem with SFRs estimated from UV data is the 
need for a extinction correction, which can be highly 
uncertain and redshift - depen d ent (flglesias-Paramo et al.l 
[2007llBurgarella et al.l [20071 ISalim et al.l I2009D . On 
the other hand, IR-based SFRs estimated by fitting 
the MIR-to-mm flux es with dust emission templates 
are model dependent (|Papovich fc Bell 2002, Dal e et all 
120051 ICaputi et al.ll2006D . Furthermore, these tracers are 
based on the assumption that the bulk of the IR emission 
traces warm dust heated by young star-forming regions. 
Thus, if a fraction of the energy heating the dust origi- 
nates from an alternative source, such as de eply dust en- 
shrouded AGNs or diff use radiation fields ()Daddi et al.l 
120071 ISalim et al.ll2009D the SFR will be overestimated. 
Nevertheless, despite these second-order effects, the un- 
certainties in the SFR are frequently driven by the ab- 
sence of sufficient IR photometry to constrain the mod- 
els robustly. In the last few years, the studies of SFRs 
at high redshift have often been based on the observed 
flux at MIPS 24 fim only and, although IR monochro- 
matic luminosities are known to correlate well to total 
IR luminosity, recent works based on more detailed IR 
coverage have demonstrated that SFRs fr om MIPS 24 im\ 
data may present significant systematics (|Papovich et al.l 
[20071 IDaddi et al.ll2007L iRiebv et alll2008D . 

In this context, iBarro et al.l 120111 (hereafter Paper I) 
presented a multi-band photometric and spectroscopic 
catalog (including data from X-ray to radio wavelengths) 
in the Extended Groth Strip (EGS), that can be used 
as a starting point for detailed analysis of the galaxy 
population. That paper describes the method used to 
measure coherent multi-band photometry and presents 
the general properties of the merged catalog, including 
an analysis of the quality and reliability of the photom- 
etry. Paper I also presents Rainbow Navigator, a pub- 
licly available web-interface that provides access to all 
the multi-band data products. 

In this paper, we focus on fitting the optical-to-NIR 
SEDs and IR emission of all the sources presented in 
Paper I using stellar population synthesis models (SPS) 
and dust emission templates. We then use the SEDs 
and fits to estimate photometric redshifts, stellar masses, 
and SFRs. We also quantify the uncertainties attending 
these estimations. In particular, we assess the quality of 
the photometric redshifts by comparing our results with 
spectroscopic redshifts and with other photometric red- 
shift compilations found in the literature. We explore 
the systematic uncertainties in the stellar masses associ- 
ated with the modeling assumptions, such as the choice 
of SPS models or the initial mass function (IMF). Fi- 
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Fig. 1. — Filter transmission for the photometric bands included 
in the dataset. The curves include the atmospheric transmission 
(for ground based observations), quantum efficiency, and the trans- 
mission of the optical elements. The curves are normalized at the 
maximum value of the transmission and scaled arbitrarily for visu- 
alization. The color code for each filter corresponds to the labels 
shown above. The optical bands depicted from top to bottom 
are those of CFHTLS, MMT and CFHT12k, respectively. The 
NIR bands, also shown from top-to-bottom are those of WIRC, 
CAHA and i/ST/NICMOS. The GALEX (FUV, NUV) and IRAC 
([3.6], [4.5], [5.8], [8.0]) filters are also listed in the top row. 

nally, we study the systematic uncertainties in the IR- 
based SFRs estimated with different IR templates and 
indicators (e.g., different total IR luminosity-to-SFR cal- 
ibrations) . 

The outline of this paper is as follows. § [2] briefly re- 
views the available data and then summarizes the most 
relevant steps of the photometric measurement and band- 
merging procedure (presented in paper I), as well as the 
overall photometric properties of the IRAC 3.6+4.5 /im- 
selected catalog. § [3] describes the techniques developed 
to perform the UV-to-IR SED fitting, and the meth- 
ods used to estimate redshifts. § [4] describes the stellar 
masses estimation technique, and quantifies the uncer- 
tainties introduced by the modeling assumptions. § [5] 
describes the methods used to fit the FIR emission to 
dust emission templates and the estimation of IR lumi- 
nosities and SFRs. § [5] presents tables containing all 
the data products presented in this paper, as well as the 
public database created to facilitate the access to these 
resources. 

Throughout this paper we use AB magnitudes. We 
adopt the cosmology Hq = 70 km _1 s~ 1 Mpc _1 , fi m = 0.3 
and Q\ = 0.7. Our default c hoice of SED modeling pa- 
rameters are: the PEGASE dFioc k, Rocca-Volmerangd 
[1991 library. alSakoeterl (119551) IMF CM i [Q. 1-1001 M^). 
and a lCalzetti et al.l (|2000) extinction law. 



2. MULTI- WAVELENGTH CATALOG 

The present work is based on the multi-wavelength 
catalog of IRAC 3.6+4.5 /im-selected galaxies in the Ex- 
tended Groth Strip (EGS; a = U h U m , S = +53°30') 
presented in Paper I. The catalog contains all the pub- 
licly available data provided by the All- Wavelength Ex- 
tended Groth Strip International Survey (AEGIS) collab- 
oration and some proprietary data including the follow- 
ing bands: GALEX FUV and NUV, CFHTLS u*g'r'i'z', 
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TABLE 1 

Photometric properties of the dataset 



Filter 




m Um [AB] 


FWHM 


Gal. ext 


Offset 


(1) 


(2) 


(3) 


( 4 ) 


(5) 


(6) 


OAT,FX-FUV 


1 Q n in 

A-iJ*J*CJ 11111 


25.6 


5.5" 


0.195 


0.04 




iiO-l.U 11111 


25.6 


5.5" 


0.101 


0.08 


MMT-m 


362.5 nm 


26.1 


1.0" 


0.049 


-0.09 


CFHTLS-?/* 

V_^l 11 _1 iJU LL 


vJO-L.-l 11111 


25.7 


0.9" 


0.045 


-0.04 


CFH-B 


^tOiy.U 1X111 


25.7 


1.2" 


0.036 


0.04 


MMT-j 


481.4 nm 


26.7 


1.3" 


0.031 


-0.09 


CFHTLS-g' 


486.3 nm 


26.5 


0.9" 


0.031 


0.03 


ACS-Vkok 


591.3 nm 


26.1 


0.2" 


0.022 


0.02 


CFHTLS-r' 


625.8 nm 


26.3 


0.8" 


0.020 


0.03 


Subaru-R 


651.8 nm 


26.1 


0.7" 


0.019 


0.00 


CFH-R 


660.1 nm 


25.3 


1.0" 


0.019 


-0.03 


CFHTLS-i' 


769.0 nm 


25.9 


0.8" 


0.015 


0.03 


MMT-i 


781.5 nm 


25.3 


1.0" 


0.015 


0.00 


ACS-igi4 


813.2 nm 


26.1 


0.2" 


0.014 


0.00 


CFH-I 


833.0 nm 


24.9 


1.1" 


0.013 


0.02 


CFHTLS-z' 


887.1 nm 


24.7 


0.8" 


0.012 


-0.02 


MMT-z 


907.0 nm 


25.3 


1.2" 


0.011 


-0.11 


NICMOS-Jno 


1.10 fim 


23.5 


0.7" 


0.008 


0.00 


tt2k - J 


1.21 nm 


22.9 


1.0" 


0.007 


-0.16 


WIRC-Jt 


1.24 /im 


21.9 


1.0" 


0.007 


0.01 


NICMOS-Hieo 


1.59 /im 


24.2 


0.8" 


0.005 


0.00 


W - K 


2.11 nm 


20.7 


1.5" 


0.003 


-0.10 


Subaru-MOIRCS-Ks 


2.15 fim 


23.7 


0.6" 


0.003 


-0.04 


WIRC-Kt 


2.16 /im 


22.9 


1.0" 


0.003 


0.00 


IRAC-36 


3.6 (im 


23.7 


2.1" 


0.001 


0.00 


IRAC-45 


4.5 /im 


23.7 


2.1" 


0.001 


0.00 


IRAC-58 


5.8 (im 


22.1 


2.2" 


0.001 


0.12 


IRAC-80 


8.0 fim 


22.1 


2.2" 


0.000 


0.12 



Note. — 

f The photometry was not measured, but taken from a published catalog. 
Col(l) Name of the observing band and instrument. 

Col(2) Effect ive wavelength of the filter calculated by convolving the Vega spec- 
trum llColina fe B ohlin 1994) with the transmission curve of the filtcr+dctcctor. 
Col(3) Limiting AB magnitude of the image estimated as the magnitude of a 
SNR— 5 detection (sec § |2. 1| for details on the flux measurement). 
Col(4) Median FWHM of the PSF in arcscconds measured in a large number of 
stars (see § 5.4 of paper I for details on the st ellarity criteria). 

Col(5) Galactic extinction estimated from the Schlcgcl ct al. (1998) maps and as- 
suming and average value of E(B-V)— 0.004. 

Col(6) Zero-point corrections applied to the photometric bands, computed by com- 
paring observed and synthetic magnitudes for spectroscopic galaxies (sec £j I3.3[t . 

GALEX, HST and MOIRCS). This region is essentially 
a field with the side edges following the contours of the 
IRAC image, i.e., inclined by 50° east of north, and upper 
and lower boundaries limited by 52.16°< 8 <53.20°. The 
bottom-right side is also restricted to a >214.04° due to 
the intersection with the CFHTLS mosaic (a square field 
oriented North up, East left). 

The 0.13deg 2 outside of the main region (hereafter 
referred to as flanking regions) also have solid optical- 
to-NIR coverage. However, the overall data quality is 
slightly lower than in the main region. The median cov- 
erage includes only 11 bands, and for the most part lacks 
the deepest, highest-resolution imaging. As a result, the 
quality of the SED coverage in the flanking regions is 
significantly lower than in the main region. For these 
reasons we focus in this contribution on the main region. 

2.1. Multi-band identification and photometry 

The procedure followed to build consistent UV-to- 
FIR SEDs from the m ultiple datasets is described de- 
tail in Paper I (see alsoJP erez-Gonzalez et"al"1l2005l and 
iPerez-Gonzalez et al.ll2008l hereafter PG05 and PG08). 
This Section summarizes the most relevant elements of 
the method, so that the impact of the photometric un- 
certainties on the parameters estimated from the SED 



MMT-u'giz, CFHT12k BRI, ACS V 6 06«8i4, Sub- 
aru R, NICMOS J lw H leo , MOIRCS K s , CAHA-JK S , 
WIRC JK, the four IRAC bands at 3.6, 4.5, 5.8, 
and 8.0 ^m, and lastly MIPS 24 and 70 ^m. We 
cross-correlated our IRAC-selected catalog with the 
X-ray (Chandra) and radio (VLA/20cm ) catalogs of 
lLaird et all (|2009l ) and llvison et al.1 (|2007h . and with all 
the spectroscopic redshifts from DEEP2 DR3 and a small 
sample of 238 spectr oscopically confirme d Lyman break 
galaxies (LBG) from lSteidel et all d20 03V The reader is 
referred to Paper I and iDavis et al.ll2007L and references 
therein, for a detailed description of all these datasets. 
Figure [Tj illustrates the different filter transmission pro- 
files for each band, and Table Q] presents the effective 
filter wavelengths, the survey depths and image quality 
achieved in each band, and the (small) zero-point re- 
calibrations (Section 13. 3p . 

The photometric coverage of the EGS is largely in- 
homogeneous, with each band covering a different por- 
tion of the IRAC mosaic (jDavis et al.ll20 07). Fortunately, 
there is a natural way to divide the field into two smaller 
sub-regions. The main region, defined by the overlap- 
ping area of the CFHTLS and IRAC frames (0.35 deg 2 ), 
presents the densest coverage (^19 bands, including 
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modeling can be assessed (§[4] and § [5j . 

First, multi-band identification is carried out by cross- 
correlating the 3.6+4.5 /mi selection with all other opti- 
cal/NIR catalogs (pre-computed with SExtractor; Bertin 
& Arnouts 1996) using a 2" search radius. The MIPS, 
Radio and X-ray catalogs required a different approach. 
For the MIPS and radio catalogs we used a 2.5" and 
3" matching radius, respectively. For the X-ray catalog 
we used a 1 or 2"radius depending on whether the X- 
ray sources were pre-identified in any other band (Laird 
et al. 2009). When two or more optical/NIR counter- 
parts separated by >l"(approximately half the FWHM 
in IRAC-3.6) are identified within the search radius, we 
apply a de-blending procedure to incorporate the mul- 
tiple sources in the catalog (e.g., irac070100 would be- 
come irac070100_l and irac070100_2). Roughly 10% of 
the IRAC sources present 2 or more counterparts in the 
ground-based images. 

Once the sources are identified, the photometry was 
computed separately in all bands, to properly account 
for the significant differences in spatial resolution. The 
fluxes were then combined to derive the merged SED. 
The procedure is carried out using our custom software 
Rainbow based on the photometric apertures obtained 
from a previous SExtractor run. 

For the optical and NIR bands, total fluxes were esti- 
mated using Kron (1980) elliptical apertures. The prop- 
erties of the aperture arc the same in all bands (although 
different between objects) and are defined from a refer- 
ence image, which is chosen by sorting the bands accord- 
ing to depth and picking the first band with a counterpart 
positive detection. Thus, this image is usually among the 
deepest, and presents a spatial resolution representative 
of the entire dataset (typically SUBARU-R or CFHTL- 
i'). Nevertheless, as a precaution, we established a min- 
imum aperture size equal to the coarsest seeing in all 
bands (1.5"). Although the choice of reference band de- 
pends on the cross-identification, the flux is measured in 
all bands independently of the counterpart detection. If 
a source is detected by IRAC only (i.e., there is no op- 
tical/NIR reference image), we use a fixed circular aper- 
ture of minimum size. If the source is detected in just a 
few optical/NIR bands (e.g., it is detected in r but not in 
z) we still use the reference aperture in the un-detected 
bands. In this way we recover fluxes for very faint sources 
not detected by SExtractor. If the forced measurements 
do not return a positive flux, the background flux from 
the sky rms within the aperture instead. These non- 
detections were not used for the subsequent SED fitting 
procedure. 

The IRAC photometry was computed using circular 
apertures of 2" radius and applying an aperture correc- 
tions estimated from empirical PSF growth curves. The 
measurement is carried our simultaneously in the four 
IRAC channels, using the 3.6+4.5 /mi positions as priors 
for the 5.8 and 8.0 /mi bands, which are much less sen- 
sitive. In the case of blended IRAC sources (i.e., those 
with multiple optical/NIR counterparts), we recomputed 
the photometry a pplying a deconvolutio n method simi- 
lar to that used in lGrazian et al.l (|2006aft or lWuvts et al.l 
(2008), which essentially relies in using smaller 0.9" ra- 
dius apertures with larger aperture corrections. Paper I 
describes the accuracy of the deblending technique. For 
the GALEX (FUV, NUV) bands we drawn the photom- 



etry from the source catalog of the public data release 
GR3. This is co mputed with aperture photometry based 
on SExtractor (jMorrissev et alj |20Q7( ). For the IRAC 
sources missed in this catalog (only ~8% and 25% of the 
IRAC catalog is detected in the FUV and NUV bands, 
respectively; see Table 4 of paper I) we used the forced 
measurement method described above. The photome- 
try in the MIPS (24 /jm, 70 /mi) bands was carried out 
using PSF fitting with IRAF-DAOPHOT and aperture 
corrections (see PG05 and PG08 for more details). 

The photometric uncertainties were computed simulta- 
neously with the flux measurement. Although the Rain- 
bow measurements are SExtractor-based, the SExtrac- 
tor photometric errors were not used, because these are 
often under estimated due to correlated signal in adja- 
cent pixels (|Labbe et al.l 120031 iGawiser et al.ll2006l) . In- 
stead, we used three different approach that range from a 
SExtract or-like method to a p rocedure similar to that de- 
scribed in Labbe et al. (2003, i.e., measuring the sky rms 
in empty photometric apertures at multiple positions). 
The photometric uncertainty was set to the largest value 
thereby derived. 

The final multi-wavelength catalog contains 76,493 
and 112,428 sources with [3.6]<23.75mag and 
[3.6]<24.75mag, respectively (these magnitude cuts 
correspond to the 85% and 75% completeness levels of 
the IRAC mosaics). Approximately 68% of the sources 
are located in the main region (52,453; [3.6] < 23.75). 
Spectroscopic redshifts have been assigned to 10% of 
the sample (only 120 are at z>1.5). A total of 2913 
stars have been identified based on several optical/NIR 
color criteria (see Section 5.4 of paper I). A source 
was identified as a star only if 3 or more criteria were 
satisfied. The stellarity value (as the total number of 
criteria satisfied) is given in Table [7] (See § [5]). The 
fractions of IRAC sources detected at 24 /im and 70 /mi 
are 20% and 2%, respectively. Finally, a total of 990 
and 590 so urces are d e tected in t he X-ray and radi o 
catalogs of lLaird et all (|2009t) and llvison et al.l (|2007t) . 
respectively. 

In the following, we analyze the SEDs and physical 
properties of the IRAC sources with [3.6]<23.75mag 
(typically SNR>10). Nevertheless, the catalog contains 
sources up to [3.6]<24.75 (3a limiting magnitude). The 
complete catalog is available in the electronic edition of 
the journal or through our web interface Rainbow navi- 
gator (see Paper I for more details). 

2.2. Galactic Extinction 

The EGS field lies at high galactic latitude benefit- 
ing from low extinction and low Galactic/zodiacal in- 
frared emission. W e derive an av e rage E(B-V)=0.004 
based on the maps of lSchlegel et al.l (|1998l ) based on sev- 
eral positions evenly spaced along the strip, centered at 
a =241.80°, S =52.80°. In our analysis, a differential 
g alactic extinction for each band is computed assuming 
a lCardelli et al.l (|1989h curve with R=3.1. These correc- 
tions, summarized in Table 1, are not included in the 
photometric catalog (presented in paper I) but these are 
applied before applying the SED fitting procedure. 

3. SED ANALYSIS: PHOTOMETRIC REDSHIFTS 
3.1. Rainbow code 
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Fig. 2. — Example of the full Spectral Energy Distribution (SED) 
of a galaxy in our sample (black dots), and fit (blue lines) of the ob- 
served UV-to-MIR photometry to a set of empirical templates com - 
puted from PEGASE 2.0 models l|Fioc fc Rocca-V olmcran g3ll997l') 
assum ing a Salpeter IMF (M S [0.1-100]M„1. and ICalzetti et ahl 
(I2000Q extinction law (see § [3j , and the FIR photometry (MIP S 
24 and 70 fim) to du st e mission models of Chary & Elbaz (2001), 
IDale fc Heloul ||2u03) and lRieke eFaLI (1200911 (see Sectionl5.lTl. The 
multiple lines in the FIR region correspond to best-fitting template 
from each of the dust emission models, and the average value of the 
three. In the upper left corner, we indicate the photometric red- 
shift, and the stellar mass, IR-based SFR and total IR luminosity 
estimated from fitting procedure. 

We computed photometric redshifts for all IRAC 
sources from the multi-color catalog presented in Paper I 
using our own dedicated template fitting code (Rainbow 
software hereafter; see PG05 and PG08). The program 
creates a grid of redshifted galaxy templates in steps of 
5z=0.01 and then applies a \ 2 minimization algorithm 
to find the template best-fitting the multi-band photom- 
etry. Upper limit detections and fluxes with uncertain- 
ties larger than 0.5 mag are not included in the fit. The 
X 2 definition takes into account the flux uncertainties of 
each band, being defined as: 



JV(band) 

\ f-^obs.- 



A- ft 



temp,! 



i=0 



(i) 



where F jj S ^ is the observed flux in the i filter and o~i 
is its uncertainty, F tem p,i is the flux of the redshifted 
template in the i filter (obtained by convolving the tem- 
plate with the filter transmission curve). A scaling fac- 
tor is applied to the input template to fit the observed 
photometry. This normalization parameter A is used to 
compute quantities such as the stellar masses, absolute 
magnitudes or SFRs (see Sections 14 . f I and [S3]) . 

Prior to the x 2 minimization procedure, the Rainbow 
code gets rid of deviant and redundant photometric data 
points. The fluxes presenting a very steep gradient with 
respect to the surrounding bands are flagged and re- 
moved before attempting the final fit. 



By analyzing the x 2 (z) distribution of the best fit in 
the model grid, we built the redshift probability distri- 
bution function (zPDF), from which we computed the 
most probable redshift and la errors, Zbcst and a z . The 
single value that minimizes x 2 ( z ) is z pcak- We found that 
Zbcst provided the most accurate results presenting less 
outliers and a smaller scatter when compared with spec- 
troscopic redshifts. The uncertainties in the photomet- 
ric redshifts are used to compute the uncertainties in the 
stellar parameters derived from the best-fitting template. 

The Rainbow code also analyzes the dust emission on 
sources with at least one flux measurement beyond rest- 
frame 8^m, i.e., the MIPS 24 and 70 ^m bands (see 
Section 15. 2p . Figure [2] shows the combined optical and 
IR SED along with the estimated physical parameters 
for a galaxy at z~l as an example of the optical and IR 
fitting techniques described here and in ^5.21 The best 
fit optical template to the data was used to estimate the 
photometric redshift, stellar mass (see Sections [3] and |4j 
and also the rest-frame UV flux. Moreover, IR luminosi- 
ties and star formation rates were obtained from the best 
fit IR template to the data at rest-frame A >5/im (see 
Section ETT]) . 

3.2. Stellar population templates 

The stellar templates used by the minimization code 
are extracted from a library of synthetic templates built 
by fitting stellar population synthesis and dust emission 
models to a representative sample of galaxies at differ- 
ent redshifts. This reference sample is drawn from the 
GOODS-N and GOODS-S IRAC surveys and have highly 
reliable spectroscopically confirmed redshifts (0<z<3) 
and at least 10 measurements of the SEDs from the UV 
to the MIR. A detailed description on the modeling of 
these templates is given in PG08. Here we briefly sum- 
marized their most relevant characteristics. 

The stellar emission of the reference template set 
was characterized using the PEGASE 2.0 models 
(|Fioc fc Rocca-Volmerangel 1 19971) assuming a Salp eter 
IMF (M e [0.1-100]M o ), and iCalzettieTal] (120001) ex- 
tinction law. We also considered the contribution from 
emission lines and the nebular continuum emitted by ion- 
ized gas. The models were obtained assuming a single 
population (1-POP models), characterized by an expo- 
nential star formation law. As a result, each template is 
characterized by 4 parameters in the 1-POP case, namely 
the time scale r, age t, metallicity Z and dust attenua- 
tion A(V). The MIR/FIR region of some templates in- 
cludes a contribution from a hot dust component that 
was computed from dust emission models using a similar 
procedure to that described in ^ 15.11 of this paper. 

Defining a representative spectral library is a critical 
issue for photometric redshift c odes, specially wh en NIR 
selected samples are studied ()Kriek et all 120081) . The 
reference sample should span a wide range of redshifts 
and galaxy colors that probe the parameter space of 
the magnitude limited sample in sufficient detail. This 
is why we included in the template set a few z>1.5 
galaxies which could not be fitted accurately with low- 
z templates. Furthermore, we complemented our syn- 
thetic templates with QSO and AGN empirical tem- 
plates drawn from iPolletta et al.l (|2007| ) that account for 
the galaxy population whose UV-to-NIR emission is not 
dominated by stars but by an AGN. 
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The template library contains a total of 1876 semi- 
empirical templates (see PG08 for more details and ex- 
amples of the SEDs) spanning a wide range of colors 
and physical parameters. Figure [3] shows that the loci of 
the observed and template colors present an overall good 
agreement for the majority of the spectroscopic galaxies 
in a wide range of optical and NIR colors. The combina- 
tion of colors based on the CFHTL S filters (panels 1 t o 
4) are consistent with Figure 2 of lllbert et al.l (|2006a), 
that presents the same colors for a sub-sample of i-band 
selected galaxies in the CFHTLS-D1 field. 

On the other hand, we find small discrepancies between 
templates and observations in the [3.6]-[4.5] IRAC color 
at low redshift (Figure [31 panel 6). This is not surprising 
considering t hat these bands are probing the rest-fra me 
NIR (see e.g. lHuang et al.ll2004tlBrodwin et al.ll2006l for 
similar examples), a wavelength region where the pre- 
dictions from stellar p opulation synthe sis models tend 
to be more uncertain (|Marastonlf2005D . Furthermore, 
these differences tend to increase at A >3 fim rest-frame, 
where galaxies can exhibit a significant contribution from 
hot dust or PAH emission features that are not con- 
templated in the optical templates and t herefore require 
more complicated modelin g procedures (Magncll i et al.l 
[20081 iMentuch et al.ll2009l) . 

3.3. Zero-point corrections and Template error function 

An improvement introduced in the current work on 
EGS over the previous Rainbow photometric redshifts 
in GOODS-N, GOODS-S and Lockman Hole (PG05, 
PG08) is the fine-tuning of the photometric zero-points 
and the use of a template error function. Both pro- 
cedures are based on the comparison of the observed 
fluxes to synthetic photometry derived from the convo- 
lution of the filter transmission curves with the best fit- 
ting templates for the galaxies wit h reliable spect r oscopi c 
redshifts. As dem onstra ted bv iBrodwin et al.l (|2006f ). 
lllbert et al.l (|2006al ) and lllbert et al.l (|2009f ). the com- 
parison between the observed apparent magnitudes and 
synthetic fluxes often shows small offsets that can lead 
to systematic errors in the calculation of the photometric 
redshifts. These offsets can be the result of small system- 
atic errors in the absolute calibration, uncertainties in 
the filter transmission curves, or they can be the result 
of intrinsic limitati ons of the templates in reproducing 
the observed SEDs ([Brammer et al.|[2008|) . 

To tackle these issues and improve the photometric 
redshift estimation, we fit the SEDs of the galaxies with 
secure spectroscopic redshift to our template set fix- 
ing the photometric redshift to the spectroscopic value. 
Then, we compute the difference between the observed 
fluxes and the template fluxes for each band, and we con- 
sider this residual value as a function of the rest-frame 
wavelength. The left panel on Figure [4] shows the result 
of applying this process to the sub-sample of galaxies 
with spectroscopic redshift and photometric fluxes with 
SNR>5 in the main region. The median of the residu- 
als (thick green line) shows an overall good agreement 
between templates and observations, with an rms (thin 
green lines) of ^2 times the median value of the photo- 
metric uncertainty across all the wavelength range (red 
lines). However, significant deviations appear in the rest- 
frame wavelengths around 200 nm, the 500-1000 nm re- 
gion and the mid-IR (A > 3 /xm). 



To diminish the effect of these discrepancies, we con- 
sidered two corrections: 1) we applied small calibration 
offsets in each band based on the residuals of the com- 
parison with synthetic magnitudes (note that these cor- 
rections refer to observed wavelengths); and 2) we used 
a template error fun ction such as that introduced in 
IBrammer et~aTl (|2008ft . 

Figure [5] shows the comparison between observed and 
synthetic magnitudes for three different i-bands (ACS, 
CFHTLS and MMT; left panel) and the u*,z',J and K 
bands (from MMT, WIRC and MOIRCS; right panel) 
as a function of redshift. The values in the parenthesis 
quote the median correction applied to each band to min- 
imize the differences with respect to the synthetic fluxes. 
Note that the three i-bands present a similar trend at 
z>l, where the observations are slightly brighter than 
the predictions from the templates. This suggest that 
the feature is related to the templates and not to the 
absolute calibration of the bands. At z>l the i-band 
(A ff ~800 nm) probes rest-frame wavelengths around 
^300-400 nm, where the overall quality of the fit to tem- 
plates is reduced. 

The overall shape of the residual distribution, shown 
in the left-panel of Figure [H is very effective for iden- 
tifying systematic deviations in the templates. This is 
because small zero-point errors in any of the individ- 
ual bands are smoothed over the rest-frame wavelength 
range due to the mixed contribution from multiple bands. 
Therefore, based on the overall scatter in the residual 
with respect to the median photometric errors, we can 
compute a template error function that parametrizes the 
overall uncertainties in the tem plates as a function o f 
wavelength. As demonstrated bv IBrammer et all (2008) , 
this function can be efficiently used as a weight term in 
the x 2 function of the SED fitting procedure to mini- 
mize the impact of the template uncertainties in some 
wavelength ranges. The bottom of the right panel of 
Figure 0] we show the median value of the absolute dif- 
ference between observed and template fluxes divided by 
the photometric error and multiplied by 0.67 to scale the 
m edian (50%) to a 1 cr(6 8%) confidence interval, as done 
in IBrammer et ai1 l2008. Compared to the results of this 
work our combination of templates and filters present a 
slightly better agreement in the rest-frame UV and NIR 
(between 1-2 /iin) , probably as a result of our larger tem- 
plate set, which present more diversity in their spectral 
shapes. 

In principle, the zero-point corrections and the effects 
of the template error function produce similar effects. 
Moreover, the re-calibration of adjacent (sometimes very 
similar) bands tend to modify the residual of both fits. 
Therefore, in order to obtain consistent results, both the 
template error function and the zero-point corrections 
are computed iteratively repeating the fitting process un- 
til we obtain variations smaller than 1-2% (typically after 
a couple of iterations). The zero-point offsets are sum- 
marized in Table [T] . Virtually all of the corrections are 
smaller than 0.1 mag, and some of them are exactly zero. 
The final results of the procedure are shown in the right 
panel of Figure [4] The application of the zero-point off- 
sets results in the flattening of the median difference be- 
tween observed and template magnitudes for the whole 
wavelength range in our SEDs except in two regions, one 
around 200 nm and the other at ~3 /im. 





rest— frame 



|>m] 



rest— frame 



[urn] 



Fig. 4. — Top: Residuals of the comparison between observed and synthetic magnitudes for a sub-sample of galaxies with spectroscopic 
redshift and SNR>5 photometry in the main region. The residuals are shifted into rest-frame wavelengths based on the effective wavelength 
of the filters and the redshift. The Figure on the left shows the raw residuals before applying the zero-point corrections nor the template 
error function to the fitting procedure. The Figure on the right shows the final result of the iterative process to compute the zero-point 
corrections and the template error function. The thick green line depicts the median value of the residuals per redshift bin. The upper and 
lower red lines indicate the median value of the photometric error at each redshift. The upper and lower thin green lines em-compasses 68% 
(lcr) of the residual distribution around the median value. Bottom: The blue line depicts the median absolute value of the residuals in the 
top panel divided by the phot ometric error and by .67 to scale the median (50%) to a 1ct(68%) confidence interval. The black line shows 
the template error function of[Branimcr ct al. (2008) divided by the median photometric error (adapted from Figure 3 of their paper). 



The poor agreement at 3 /im is most probably associ- 
ated with limitations in the NIR-MIR range of stellar 
population templates and the contribution from PAH 
emission, which is also not properly taken into ac- 
count in the SPSs models nor the dust emission tem- 
plates. In addition, there is a small peak/bump at 
^350-450 Dm, which i s very sim ilar to the feature re - 
ported bv iWild et ail ([2001 and lWalcher"et~aT] (|2008h . 
In these papers, they explain this effect with an ex- 
cess in the strength of th e B aim er break in the mod- 
els by iBruzual fc Chariot! (|2003l ). relative to the ob- 
served values. The pe ak at 200 nm c an be partly re- 
lated with the use of a ICalzetti et al.l (|2000[ ) extinction 
curve in the modeling of the the galaxy templates. This 
parametrization lacks the silicate absorption at 2175 A , 
which appears in oth er extinction curves, such as that of 
ICardelli et al.l (|1989f ) , which has been claimed to be pro- 
duced by PAHs. The presence of this absorption bump 
has been reported on some studies of high-z galaxies 
(|Noll et al.ll2007t iNoterdaeme et al.ll2009D . 

As an additional check of the accuracy of the method 
we compare our photometry against the fluxes of a con- 
trol sample of ^300 bright unsaturated stars in common 
with the SDSS. In particular, we restrict the compar- 
ison to relatively blue sources (it'-g=1.2, in the MMT 
bands) in order to avoid large color corrections in the 
filter transformations. These color terms were computed 
by convolving the filter t ransmissions w ith the spectra of 
F, G and K class stars (|Kuruczlll992f ). which makes up 
for most of our sample of stars. The transformation with 
respect to the MMT bands, which present a filter system 
similar to that of SDSS, are 



'MMT — *SDSS ' 
ZMMT = ZSDSS 



0.203 • [i - 
- 0.087- \i 



*]sdss - 0.002 
■ *]sdss - 0.002 



After applying these corrections, we find zero-point 
offsets of Au'=-0.05,Ag=-0.10,Ai=-0.01 and Az=-0.09 
with respect to the SDSS. The values are roughly consis- 
tent with our previous results based on galaxy templates. 
Only the u' and z bands present slightly lower values of 
the correction. These could be an additional effect of the 
template uncertainties (at 250 and 450 nm rest-frame), 
specially for the u' band .Also, it is worth noting that 
the zero-point offsets are estimated simultaneously and 
iteratively for all bands whereas the comparison to SDSS 
is done separately for each band. 

3.4. Photometric redshift accuracy 

In this section we analyze the overall accuracy of the 
photometric redshifts (.Zphot) by comparing them against 
spectroscopic redshifts (z S pec)- I n particular, we study 
the quality of our results as a function of the spectro- 
scopic redshift and the observed magnitude in optical 
and NIR bands, and we provide specific results for dif- 
ferent groups of galaxies such us X-ray, MIPS or Radio 
sources. For the 76,936 galaxies ([3.6]< 23.75) in the 
sample we identify 7,636 (^10%) spectroscopically con- 
firmed sources from the DEEP2 catalog (mostly at z<l) 
and from a s mall sample of LBGs (z~3) presented in 
iSteidel et al.l (|200l . 



3.4.1. 



DEEP2 sample 
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MMMT = w S dss - 0.095 • [li - c?]sdss 
ffMMT = fl'SDSS — 0.063 ■ [g — i]sdss 
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2 P hot versus z s] 

Figure [S] shows the comparison between z p hot and z s - 
for 6,191 and 1,445 sources with reliable spectroscopic 
redshift in t he main and flankin g regions, respectively. 

Following lllbert et al.1 (|2006bl ) . we quantified the red- 
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Fig. 5. — Left - Differences between observed and synthetic magnitudes as a function of redshift in the ACS-«814, CFHTLS-i' and MMT-i 
bands. The values quoted in the parenthesis indicate the zero-point correction applied to these bands. The dashed blue line depicts the 
median difference between observed and synthetic photometry after the zero-point correction has been applied. The green lines show the 
median photometric uncertainty in each band as a function of redshift multiplied by a factor 2. At z~l the residuals are dominated by 
a systematic offset in the templates instead of by deviations in the photometric calibration. . Right - Same as the left panel for the 
MMT-u*,MMT-z',WIRC-J and MOIRCS-K. bands. 



shift accuracy using the normalized median absolute de- 
viation (<7 N mAd) Of Az = Zp ho t - z s 



~spcc 



cnmad = 1-48 x median 



Az — median(Az) 



1 + z s 



(2) 



This quantity is equal to the rms for a Gaussian dis- 
tribution and it is less sensitiv e to the out l iers th an the 
usual rms divided by (1+z) (jllbert et al.ll2006bO . We 
define r\ as the fraction of catastrophic outliers (those 
sources having |Az|/(l + z) >0.20). 

Table [2] and [3] summarize the quality of z p hot as a 
function of redshift in the main and flanking regions. 
The overall scatter and median systematic deviation 
are ct N mad=0.034 and 0.046, and A«/(l+z)=0.010 and 
0.013 for each region, respectively. As expected, the rms 
in the flanking regions, where the overall photometric 
quality is slightly lower, is higher (~20%) than in the 
main region. Nonetheless, the outlier fraction is only 1% 
worse. 

The bottom panels of Figure [5] show the density plot 
of Az/(1 + z) as a function of redshift. The subset of 
LBGs at z spcc >2.5 are shown as dots. These sources are 
explicitly discussed in the following section. The scat- 
ter distribution indicates that the accuracy of z p hot does 
not depend strongly on the redshift up to the limit of 
the DEEP2 sample. The systematics in both regions are 
fairly similar presenting a minimum scatter at 0.5<z<l, 
around the peak of the z p hot distribution (see § 13 - 7[) , and 
increasing by a factor of ~1.3 at lower and higher red- 
shifts (zspcc <l-5). We find that the slightly worse perfor- 
mance at z<0.5 is associated to the use of 4 IRAC bands 
in the fitting of the SEDs. Although the template error 
out-weights the contribution of these bands (mostly at 
A >3 fjxa rest- frame; see § I3.3[) their contribution cause a 
broadening of the zPDF that tends to increase the scat- 



ter. Nonetheless, this effect does not increase the outlier 
fraction at z<0.5, which is comparatively lower than at 
l<z<1.5, for similar values of ctnmad- 

We also analyze the quality of z p hot as a function of 
the optical and NIR magnitudes. As the efficiency of 
•z p hot mostly relies on the detection of strong continuum 
features, the estimates are highly sensitive to overall con- 
sistency of the multi-band coverage. Figure [7] shows the 
scatter in Az/(l+z) as a function the observed magni- 
tudes in the R and [3.6] bands for sources in the main 
region. The results in the flanking fields are similar, but 
with a larger scatter. We choose these bands to be rep- 
resentative of the brightness of the sources in the op- 
tical and NIR, and ultimately of the overall band cov- 
erage. Note that, although this is NIR selected sam- 
ple, most of the photometric coverage consist on optical 
bands. Thus, galaxies with faint optical magnitudes tend 
to present worse photometric redshifts. The magenta 
bars depict the median deviation and ctncmad per mag- 
nitude bin. We have corrected both plots by a median 
offset of Az/(l+z)=0.01. In the R-band, the scatter in- 
creases monotonically with the optical magnitude from 
(Tncmad = 0.03 - 0.06 for R=22-25, and >50% of the 
outliers are located at R>23.5. The scatter is also wider 
at R<22. However, since most of these bright galaxies lie 
at low-z, this trend is essentially the same one mentioned 
above for sources at z <0.5. Interestingly, there is weaker 
dependence in the scatter (and the outlier fraction) with 
the [3.6] magnitude than with the R-band magnitude. 
This is because, the overall quality of the optical pho- 
tometry is more relevant for constraining the shape of 
the SED, and there is typically a wide range of optical 
brightnesses for any given [3.6] magnitude (see, e.g., Fig- 
ure 6 of Paper I) . 



3.4.2. 



2-phot 



versus z sp cc-' LBGs sample 



10 




N 
+ 



N 




N 



4 
3.5 

3 
2.5 

o 

0) 

9- 2 
1.5 

1 

0.5 



FLANKING REGIONS 
1 445 galaxies 
[3.6]<23.75 
r|=3% 

o A ,,, .= 0.046 

Az/(1+z) 

PLGs 
X-ray 



N 
+ 



N 

< 



0.5 



1.5 2 
^phot 



2.5 



3.5 




Fig. 6. — Top panels: Photometric redshifts versus spectroscopic redshifts for [3 .6l<23.75 m a g sou rces in the main region (left) and 
flanking regions (right). Blue stars indicate sources detected in the X-ray catalog of Laird ct al. (2009). Red dots correspond to galaxies 



with a power-law spectrum in the IRAC bands. This feature is frequently used to identify obscured AGNs ( Alonso-Hcrrcro ct al. 2004, 
IDonlev ct al. 2007), that usually under-perform in the photometric redshift procedure. Bottom plots: Density plots of the scatter in Az/ (l+z) 
as a function of redshift for the main region (left) and flanking regions (right). Each contour contains (from the inside out) 25%, 50%, 75% 
and 90% of the spectroscopic sample, respectively. 



Given that the DEEP2 spectroscopic catalog consist 
mostly of low redshift galaxies (68% is located at z<0.9), 
we have included in our sample spectroscopic redsh ifts 
drawn from the LBG catalog of ISteidel et al.l (|2003h to 
specifically study the accuracy of 2 p hot beyond the clas- 
sic spectroscopic limit. This catalog contains 334 LBGs 
galaxies, 193 of them with confirmed spectroscopic red- 
shift. To check the quality of our z p hot at z>2.5, we first 
compare our results to their z speC j and then we check that 
the our photometric redshift distribution for the whole 
LBG sample is consistent with the average redshift of 
this population. 

We identify IRAC counterparts for 91(147) of the spec- 
troscopic LBGs with [3.6]<23.75mag (24.75). The rest 
were missed mainly because they lie out of the observed 
area in the IRAC survey; only 10 galaxies were lost 
due to their faintness in the IRAC bands. Note that, 
although these LBGs are relatively bright in the opti- 
cal (i?<25.5mag), most of them are intrinsically faint 
in the IRAC bands, ^50% and 20% are fainter than 
[3.6]=23.75mag and 24.75 mag, respectively. In gen- 
eral, LBGs are known to span a wide range of IRAC 
magnitudes (jHuang et al.ll2005llRigopoulou et alJ l2006). 
and they exhibit a clear dichotomy in the R-[3.6] color, 
with red (R— [3.6] >1. 5) source s showing brighter IRAC 
magnitudes than blue sources (|Magdis et al.ll2008h . We 
find that the median magnitudes and colors for the 
LBGs in our sample are [3.6]=22.74, R-[3.6]=2.06 
and [3.6]=23.80, R-[3.6]=0.88 for red and blue galax- 



ies, respectiv el y, in good agreement with the values of 
iMagdis et ail (|2008l ) for a larg e sample of LBGs also 
drawn from the LBG catalog of ISteidel et al.l (|2003l ). 

The quality of £ph t f°r the spectroscopic LBGs is sum- 
marized in Table [2] For the galaxies with [3.6] < 2 3. 75, 
both the scatter and the outlier fraction (oncmad = 
0.063, rj = 10%) are slightly worse than the median of 
the sample, as expected by their intrinsic faintness in sev- 
eral optical and NIR bands. Nonetheless, the statistics 
are similar to the results of o ther authors at high red- 
shift (e.g., IWuvts et al.ll2008l ) indicating that our z pno t 
still provide reasonably consistent values beyond z>1.5. 
If we also consider the faintest sources ( [3. 6] < 24. 75), the 
statistics do not degrade much (ctncmad = 0.069), even 
though we are including 60% more sources. We have 
visually inspected the outliers and at least 4 of them 
present flux contamination from close-by sources and an- 
other 3 are strong AGNs detected in the X-rays. The 
rest of them present a high-z solution in the zPDF, but 
the flux at [5.8] and [8.0] is too faint to reliably identify 
the rapid decline of the stellar component at A >1.6/im, 
which results on favoring the low-redshift solution. 

We also compared the photometric redshift distribu- 
tion of the 155 galaxies with [3.6]<23.75 identified in the 
whole LBG catalog. Figure [5] shows the redshift distribu- 
tion of the photometric and spectroscopic LBGs in our 
sample. The median value and quartiles for the pho- 
tometric LBGs with [3.6]<23.75 is z p hot=2.8±Q g consis- 
tent with the median redshift of the spectroscopic sample 
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(2=2.95) and with the typical width of th e redshift dis- 
tribution for the LBG c riteria (z=3.0±0.3: ISteidel et al.1 
12004 iReddv et alj|2005l ). About 14% of these sources lie 
at redshift z<1.5, similarly to the outlier fraction of the 
spectroscopic sample. 



3.4.3. 



'phot 



versus z 



ipcc : X-ray, Power-law, MIPS and 
Radio galaxies 



Fig. 7. — Density plot of the scatter in Az/(1 + z) as a function of magnitude in the R (left) and [3.6] (right) bands. The magenta bars 
depict the median value of Az/(1 + z) and ctncmad (with respect zero) for each magnitude bin. The lower panel of each plot shows the 
fraction of photometric redshift outliers (rj) as a function of magnitude. 

outlier fraction than the X-ray sources as a consequence 
of having their SED more contaminated by the AGN 
emission. 

Similarly to the PLG, wc find that sources with red 
colors in the IRAC bands (f[3.6]<f[4.5]<%.8]<f[8.o])j but 
not strictly satisfying the PLG criteria, make up for up 
to 15% of the total outliers. This is not surprising given 
that, for the typical galaxy at z<2, the presence of the 
stellar bump (at 1.6 /im) causes the flux in the last two 
IRAC bands to be lower than in the previous two. Thus, 
for these sources, the code will try to assign incorrect 
high redshift values of z p hot- 

Table [5] and [3] also quote numbers and z p hot statis- 
tics for the spectroscopic sub-samples of galaxies de- 
tected in MIPS 24 ^m (f(24)>60 /iJy) , MIPS 70 fj,m 
(f(70)>3.5 mJy) a nd in the catalog of Radio sources of 
llvison et al.l J2007). The latter present a slightly worse 
accuracy than the median of the sample, whereas the 
MIPS detected galaxies present essentially the same qual- 
ity as the rest of the spectroscopic galaxies. This indi- 
cates that for most of them the IR emission does not con- 
tribute significantly to the NIR-MIR region fitted with 
the optical templates. 



We analyze in detail the quality of the z p hot for sam- 
ples of galaxies that are known to present particularly 
different SEDs from the majority of the templates (e.g., 
X-ray or AGNs), which could cause a degradation of the 
redshift estimate. These sources are shown with different 
makers and colors in Figure [6] and their z p hot statistics 
are summarized in Table [2] and 03 

The blue stars show galaxies identified in the X-Ray 
Chandra/ ACIS catalog in EGS (jLaird et al.ll2009D . prob- 
ably indicating the presence of an AGN. The SED of 
these sources is likely affected by the AGN emission, 
which in principle should decrease the efficiency of the 
template fitting procedure. In spite of showing a larger 
outlier fraction (particularly at z>1.5), the z p hot for X- 
ray sources are quite accurate, with a scatter similar to 
that of the full sample. 

The red dots in Figure [6] depict galaxies satisfying the 
power-law crit eria (PLG) commonly used t o identify ob- 
scure d AGNs (jAlonso-Herrero et al]|2004 iDonlev et all 
2007;. a good fraction of them being undetected in the 
X-rays. We find a surface density of 0.26arcmin~ 2 for 
PLGs, in good agree ment with the 0.22 arcmin -2 given in 
IDonlev et al.l ({2007, wc apply a similar criteria restricted 
to P x >0.1 and a slope a <-0.5). However, less than 2% 
of these sources present a spectroscopic redshift. Com- 
paratively, PLGs present a lower accuracy and higher 



3.4.4. Error analysis 

The lcr uncertainty of the photometric redshifts, 
Az p hotj is computed from the zPDF as the semi- width of 
the redshift range corresponding to a 68% confidence in- 
terval around the probability peak. This value allows to 
provide an estimate of the accuracy for sources without 
a spectroscopic redshift, which are >90% of the sample. 

Table [2] and [3] quote the values of Az p hot as a func- 
tion of redshift in the main and flanking regions. Based 
on these results, we find that 62% (approximately la) of 
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TABLE 2 

Photometric redshift accuracy in the Main region [3.6]<23.75 



Rainbow EAZY 



Redshift 


No. 


Az/(l+z) 


^NCMAD 






Az/(l+z) 


CT NCMAD 


ij 




Qz 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


All 


6191 


0.010 


0.034 


2% 


0.068 


0.019 


0.029 


3% 


0.059 


94% 


0.<z<0.5 


1637 


0.015 


0.040 


2% 


0.070 


0.019 


0.031 


1% 


0.049 


97% 


0.5<z<1.0 


3171 


0.007 


0.028 


2% 


0.061 


0.018 


0.025 


2% 


0.055 


95% 


1.0<z<2.5 


1292 


0.017 


0.035 


5% 


0.083 


0.021 


0.032 


4% 


0.078 


89% 


z>2.5 (LBGs) 


91 


-0.023 


0.063 


10% 


0.110 


-0.014 


0.043 


15% 


0.107 


42% 


z>2.5 (LBGs [3.6]<24.75) 


147 


-0.027 


0.069 


12% 


0.105 


0.012 


0.060 


23% 


0.119 


33% 


X-ray 


142 


0.003 


0.038 


10% 


0.070 


0.008 


0.034 


10% 


0.057 


82% 


PLCs 


8 


0.031 


0.142 


50% 


0.092 


0.018 


0.094 


25% 


0.108 


50% 


MIPS-24^m 


1955 


0.010 


0.033 


3% 


0.068 


0.023 


0.026 


3% 


0.055 


94% 


MIPS-70/im 


262 


0.015 


0.045 


1% 


0.071 


0.031 


0.028 


2% 


0.050 


95% 


Radio 


85 


-0.001 


0.052 


5% 


0.066 


0.017 


0.032 


2% 


0.048 


92% 



Note. — Photometric redshift quality with the estimates with Rainbow and EAZY (see ij I3.5.2H . 

(1) Spectroscopic redshift range. 

(2) Number of sources in the redshift bin. 

(3.7) Median systematic deviation in Az/(l+z); Az = z p hot ~ z ape c. 

(4.8) Normalized median absolute deviation. 

(5.9) Percentage of catastrophic outliers (|Az|/(l + z) >0.20). 

(6.10) 68% confidence interval in the zPDF around the most probable z p ^ ot . 

(11) Percentage of the sources with Q 2 >1 in EAZY (high quality flag according to Brammer et al. 2008). 



TABLE 3 

Photometric redshift quality in the Flanking regions [3.6]<23.75 



Rainbow 



EAZY 



Redshift 


No. 


Az/(l+z) 


CNCMAD 


') 


Az phot 


Az/(l+z) 


°"NCMAD 


V 


Az phot 


Qz 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


All 


1445 


0.013 


0.046 


3% 


0.065 


0.027 


0.050 


4% 


0.073 


87% 


0.<z<0.5 


373 


0.021 


0.058 


2% 


0.065 


0.037 


0.058 


5% 


0.052 


83% 


0.5<z<1.0 


785 


0.015 


0.040 


3% 


0.059 


0.014 


0.069 


8% 


0.077 


75% 


1.0<z<2.5 


274 


-0.001 


0.058 


8% 


0.082 


0.009 


0.043 


8% 


0.079 


70% 


X-ray 


33 


-0.019 


0.031 


9% 


0.057 


0.000 


0.035 


9% 


0.061 


88% 


PLCs 























MIPS-24^m 


416 


0.011 


0.046 


4% 


0.061 


0.028 


0.048 


5% 


0.068 


88% 


MIPS-70Atm 


66 


0.008 


0.050 


2% 


0.063 


0.026 


0.052 


3% 


0.055 


89% 


Radio 


13 


0.000 


0.055 


8% 


0.071 


0.035 


0.015 


0% 


0.060 


85% 



Note. 



Same as Table [2] 



The me- 



ii nofiO.0.56 
U - UOD 0.021 



the galaxies present values of Az p hot < |A2 
dian value and quartiles of Az p hot/(l + z) = 
in the main region are consistent with the statistics for 
oncmad and also with \Az\/(l + z) = 0.027g;gf §. A sim- 
ilar agreement is found for the sources in the flanking 
regions. Note that, as the Az p hot is computed from the 
zPDF its minimum value is limited by the step size of 
the redshift grid (Az=0.01), and thus it tends to present 
larger values than |Az|, specially for very accurate 2 p i,ot- 
Therefore, it is not surprising that Az p hot/(l + z) is on 
average larger than all the other scatter estimates. In 
fact, this indicates that Az p hot provides a robust esti- 
mate of the uncertainty in z p h t, which can be underes- 
timated if it assumed t o be equal to ctncmad (see e.g., 
ICardamone et al.ll20ldh . 

In order to obtain a better characterization of the 
catastrophic outliers caused by a poor fit to the data, 
we analyze the distribution of sources as a function of 
the reduced \ 2 °f the SED fitting. Figure [9] shows the 



distribution of —log(x 2 ) for the full photometric sam- 
ple, the spectroscopic sample and the catastrophic out- 
liers. Approximately 83% and 94% of the galaxies in 
the photometric and spectroscopic sample present values 
of x 2 lower than the median, of the outlier distribution 
(—log(x 2 ) < ~0.6), i.e., half of the outliers are located 
within the ^20% and 5% of the sources in the photo- 
metric and spectroscopic samples with the worse values 
of x 2 - 

Finally, we also find that 58% of the sources with 
significantly different values of Zbest and z pC ak (l^best- 
Zpeak|/(l+z)>0.2) are outliers. These sources account 
for only 1% of the spectroscopic sample, but they rep- 
resent ~12% of the outliers. Therefore, the difference 
between both values is another useful indicator of possi- 
ble outliers. 

3.5. Comparison of photometric redshift catalogs 

Here we compare the z p hot computed with Rainbow 
to other previously published z p hot catalog and to the 
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redshift 

Fig. 8. — Photometric redshift distribution of the photometric 
(155, red dashed line) and spectroscopically confirmed (91, black 
line) LBG s (pLBG.zLBG) wi th [3.6]<23.75 in common with the 
sample of Steidel et al. (2003|). The spectroscopic redshift distri- 
bution for the zLBG is shown as a filled grey histogram for com- 
parison. 
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Fig. 9. — Distribution of the normalized x 2 resulting from the fit 
of the data to the templates during the calculation of 2 p hot- The 
black line is for the full sample ([3.6]<23.75), and the magenta line 
and grey area are for the spectroscopic sample (x5) and photomet- 
ric redshift outliers (»7X50), respectively. The corresponding lines 
depict the median value and quartiles of each distribution. 

estimates obtained with a different code. The alternative 
Zphot are also included in the our data release (see § [S]) 
in Table 0. 

3.5.1. Rainbow vs. Ilbert et al. (2006a) 

We compare the z D hot prese n ted in this paper with 
those derived by lllbert et"aTI (|2006al hereafter, 106) 
based on optical data from the CFHTLS. These authors 
used a i'-band selected sample with i' <24 and obtained 
photometric redshifts for the four CFHTLS deep fields. 
The Zphot were computed using t he template fittin g code 
Le phare (Arnouts & Ilbert; e.g.. lllbert et al.| [2009') for ~ 



500, 000 sources observed in 5 bands u*,g',r',i',z'. Their 
template librar y is based on an upgr ade of the empiri- 
cal tem plates of IColeman et al] (|1980t ) and iKinnev et al] 
(1996) computed by applying zero-point corrections and 
interpolating between spectral types. Their z p hot also in- 
clude a Bayesian prior on the redshift distribution. The 
accuracy of their results for the D3 field (the EGS) is 
cncmad = 0.035 with r] = 4% for sources wi th i' <24 
and z<1.5. More recently. ICoupon et al.l (|2009|) repeated 
essentially the same exercise using the latest data release 
of the CFHTLS T004, obtaining z p hot of almost identical 
quality 

The source density in the i'-band selected sample of 
106 is 25, 42 and 96 sources/arcmin 2 up to limiting mag- 
nitudes of i' =24, 25 and 26.5 (the estimated SNR~5 
level). The source density of the IRAC selected catalog 
is ^44 sources/arcmin 2 at [3.6]<23.75. This means that 
their source density at i' <25, which is essentially the 
spectroscopic limit (R=25), is similar to ours. However, 
at i' <24, the limiting magnitude for their best perform- 
ing Zphot, the source density in i' is approximately 50% 
that in IRAC. At the faintest optical magnitudes, the 
source density in the z'-band selected catalog is larger, 
although the quality of these z p hot is worse than for the 
i'<2A sample, given that many of the galaxies will also 
be undetected in the shallowest optical bands (u*,z'). 

Even presenting similar source densities, the nature of 
the galaxies in an z'-band and an IRAC-selected sam- 
ples is different, and some of the sources in one selection 
will be missed by the other. We find that the optically 
bright galaxies missed by the IRAC catalog ([3.6]>23.75) 
present a median and quartile redshifts z p hot — 1-Oo'g, 
while the infrared bright galaxies undetected in the op- 
tical (?'>26.5) present z p hot =1-8|;3- The high-z sources 
missed in the IR selection are typically low-mass galaxies 
(similar to LBGs), i.e., our catalog favors the detection 
of high-z massive galaxies, as expected. 

We cross-correlated the catalog of 106 to the IRAC se- 
lected sample using a search radius of 1.5". Due to small 
differences in the extraction of the catalogs, the compar- 
ison is restricted to a slightly smaller portion of the main 
region (214.09°<a<215.72° and 52.20°< 6 <53.16°). 
Out of the 49605 IRAC sources, 40% and 88% are de- 
tected in 106 to i' <24 and 26.5, respectively. The cross- 
match to the DEEP2 spectroscopic redshifts contains 
5454 galaxies simultaneously identified in all three cat- 
alogs ([3.6]<23.75, i'<26.5). Approximately 6% of our 
spectroscopic subsample is missed due to a more conser- 
vative source removal around bright stars in 106. 

The top panels of Figure [10] show the comparison of 
z p hot versus z spec for the galaxies in common between the 
Rainbow, 106 and DEEP2 catalogs with [3.6]<23.75 and 
z'<25, without any other requirement of band coverage. 
Table [4] summarizes accuracy of the z p hot hi Rainbow and 
106 for these sources as a function of redshift. We also 
list the fraction of catastrophic outliers in each catalog 
that is recovered in the other (shown as green dots in 
Figure [10]). 

The overall scatter in 106 for the sources in common 
with the IRAC sample is consistent with their results for 
the whole D3 sample. The comparison as a function of 
redshift indicates that 106 estimates at lower redshift are 
slightly more precise, probably as a result of the template 
optimization algorithm and the Bayesian prior (see Fig. 6 
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TABLE 4 

Rainbow z p h t versus 106 at i'<25 and [3.6]<23.75 



Rainbow 106 



Redshift No. ctncmad V R(n) ctncmad V Rfa) 
0) (2) (3) (4) (5) (6) (7) (8) 

All 5454 0.034 2% 82% 0.036 5% 55% 



0.<z<0.5 1444 0.040 2% 83% 0.032 5% 43% 

0.5<z<1.0 2787 0.028 2% 80% 0.031 3% 53% 

1.0<z<2.5 1143 0.035 4% 80% 0.054 8% 63% 

z>2.5 (LBGs) 80 0.063 9% 91% 0.345 46% 42% 

Note. — 

Photometric redshift quality in the estimates with Rainbow and in 106. 

(1) Spectroscopic redshift range. 

(2) Number of sources in the redshift bin. 

(3.6) Normalized median absolute deviation. 

(4.7) Percentage of catastrophic outliers (|Az|/(l + z) >0.20). 

(5.8) Fraction of catastrophicc outliers in the other code presenting an accu- 
rate z phot . 



of 106), but also because of our slightly lower performance 
at z<0.5. On the contrary, the fraction of catastrophic 
outliers in 106 is larger than in Rainbow for all redshifts, 
and particularly at z> 2.5. Moreover, the z p hot Rain- 
bow is able to recover ^80% of these outliers. At z>l, 
our larger band coverage, mostly in NIR bands, provides 
more accurate estimates. Note that the IRAC fluxes for 
the LBGs play a critical role on providing more accurate 
redshifts (and stellar parameters) for these sources. 

3.5.2. Rainbow vs. EAZY 

Here we check again the quality and overall consistency 
of our SEDs and z p hot by computing an independent es- 
timation of the Zphot with a different code. A successful 
result using a different fitting code based on different 
template sets would certify that the catalog reproduces 
accurately the observed SEDs and is therefore suitable 
for galaxy population studies. 

We computed alter native z D hot using the photometric 
redshift code EAZY (|Brammer et al.ll2008[) . The advan- 
tage of EAZY is that it was conceived to provide ac- 
curate photometric redshift estimates for NIR selected 
samples in absence of a representative calibration sam- 
ple of spectroscopic redshifts. The code makes use of a 
new set of templates computed from a if-limited subsam- 
ple of the Millennium S imulation (jSpringel et al.l [20051 
iDe Lucia~fc BlaizotJ 120071) and modeled by fitting the 
synthetic SEDs with PEGASE models and applying an 
optimization algorithm. The final result is set of 6 tem- 
plates that essentially reproduces the principal compo- 
nents of the catalog. Furthermore, a template error func- 
tion was introduced to account for systematic differences 
between the observed photometry and the template pho- 
tometry at different wavelengths. After trying different 
configurations for the input parameters, we find that the 
best results in the z p hot-z S pcc comparison are obtained 
using the template error function and incorporating a 
B ayesian p r ior on the redshift distribution s imilar to that 
of lBemtezI ()2000L see lBrammer et al.1 [20081 for more de- 
tails). The use of the template error function is decisive 
to weight the contribution of the IRAC bands at lower 
redshifts as we have also verified in our own z p h t (see 

§& 

In addition, we find that a critical issue to avoid a 



severe contamination from catastrophic z p h t is the use 
of the purged photometric catalog produced by Rainbow. 
Prior to the fitting procedure, Rainbow carries out a first 
pass on the catalog where potential photometric outliers 
are removed. If we use the resulting catalog as input 
for EAZY the outlier fraction is reduced by a factor ^5, 
illustrating the relevance on the photometric errors not 
only in the overall quality of the z p h t but also in the 
catastrophic errors. 

The bottom panel of Figure [10] depict the comparison 
the Zphot with EAZY versus spectroscopic redshifts for 
the galaxies in the main region (i.e., the same galaxies 
depicted in the left panel of Figure [6]). As in the previous 
section green markers indicate outliers in z p h t Rainbow. 
Table [2] and [3] summarizes the quality of z p hot for Rain- 
bow and EAZY in different redshift bins for the main 
and flanking regions, respectively. The overall scatter 
and outlier fraction of both estimates are roughly similar, 
with the estimates of EAZY presenting a slightly higher 
accuracy in the main region, but lower in the flanking 
fields. We also note that z p hot with EAZY perform bet- 
ter at z<0.5, being less sensitive to the mild broadening 
of the zPDF present in Rainbow. The tables also quote 
the 68% confidence intervals of z p hot, which are similar 
for both codes and are in go od agreement with othe r 
results based on EAZY (e.g., ICardamone et al.ll2010j) . 
We show that the fraction of sources with a quality pa- 
rameter Q Z <1 in E AZY (good photometric redshifts; 
iBrammer et al.l [20081 ) is typically >90% except for the 
highest redshift bin, where the photometry is more un- 
certain, due to the intrinsic faintness of these sources, 
and thus there is larger fraction of outliers. In addition, 
we find that, roughly 40-50% of the outliers in EAZY 
present a poor z p hot estimate in Rainbow and viceversa. 
Therefore, similarly to the galaxies with different values 
of z pca k and z pro b, galaxies with significantly different 
estimates both catalogs are frequently (^50%) outliers. 

From the good agreement between the different z p hot 
estimates we conclude that the photometric catalog pro- 
vides accurate SEDs suitable for studies of galaxy popu- 
lations irrespective of the code used for the analysis. 

3.6. Number densities and redshift distribution of 
NIR-selected galaxies. 
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Fig. 10. — Top: Comparison of z p hot versus z Bpcc for the estimations presented in this paper (left) and the ones in 106 (right). The 
sample is drawn from the overlap region between the CFHTLS-D3 area and the IRAC mosaic (main region). Both figures contain the same 
5454 sources simultaneously detected in both catalogs and in the DEEP2 sample at [3.6]<23.75 and i'<24.5. Green points depict z p i lot in 
Rainbow for the outliers in 106 (left) and viceversa (right). Bottom: Comparison of z p hot computed with EAZY versus z sp cc for sources 
in the main region. This Figure is equivalent to the left panel of Figure \E\ for estimates based in EAZY. Green dots depict catastrophic 
outliers in the z p h t estimated with Rainbow. 

As an additional test of the accuracy of z p hot we com- 
pare the number densities and redshift distributions of 
NIR color-selected populations with the results from 
other authors. In order to facilitate the comparison to 
the references, the magnitudes in this section are given 
in Vega system. 

Given the highly non-uniform band coverage of the 
field, we have chosen to compute galaxy colors based on 
synthetic magnitudes. An advantage of this method is 
that synthetic photometry behaves better than directly 
observed values when deep data is not available in some 
of the required bands, allowing us to assign robust fluxes 
for undetected sources in the shallower bands. This is the 
same procedure that we used in PG08, and is similar to 



that p resented in lGrazian et al.l (|2007t ) and lQuadri et al.l 
(2007). We restrict the analysis in this section to the 
0.35 deg 2 of the main region which count with better pho- 
tometry. 

For obvious reasons, the success of this method de- 
pends critically on the quality of the synthetic fluxes. In 
§ 13.31 we showed that these fluxes provide an accurate 
representation of the observed values in the magnitude 
range covered by the observations. The median offsets 
are very small and the scatter is consistent within a factor 
~2 with the photometric errors at different magnitudes 
(see Figure [5]) . 

In order to avoid possible selection effects, we restrict 
the analysis to NIR selected galaxies which would be 
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Fig. 11. — Galaxy number counts of our sample in the _K"-band, derived from synthetic observed magnitudes, compared to results from 
the literature. The black dots connected with a solid line depict the values derived in the present work. The other symbols show the results 
from other authors. Top-left; Number counts for the complete galaxy sample ([3.6]<23.75) in the main region. Top-right; Number counts 
for DRG. Bottom-left; Number counts for s-BzK galaxies. Bottom-right; Number counts for p-BzK galaxies. 



fully represented in the IRAC-select ed sample. We se- 
lected Distant Red Galaxies (DRG; iFranx et all 120031 ) 
as galaxies with [J—K] >2.3, and BzK galaxies, both 
star forming (s-BzK) an d passively evolvin g (p-BzK), 
following the equations in lDaddi et a l. (2004). Both cri- 
teria were proposed to target massive galaxies at z— 2, 
although DRG and p-BzK are best at selecting galaxies 
with a significant fraction of evolved stars, whereas s-BzK 
select star-forming galaxies similar to those found by the 
low-redshift equiv alent of the LBG criteria (LBG/BX; 
ISteidel et al"1l2004l ). For the DRG we convolved the tem- 
plates with the VLT/ISAAC J and K filters, whereas 
for the BzK we used VLT/FORS B, HST/ACS z and 
VLT/ISAAC K w hich are the same filters used in 
iDaddi et alJ (|200l . 

The top left panel of Figure [TTJ shows the AT-band 1 
number counts for the IRAC-selected catalog compared 
to other results from the bibliography. O ur counts are 
in ver y good agreement with the values of iQuadri et al.l 
(2007) for the MUSYC survey, and with our previous 



results in Barro et al. (2009) for the South region of 
EGS (-30% overlap with the IRAC sample). The overall 
agre ement with the count s of the Palomar-WIRC cata- 
log (jConselice et all 12008ft is slightly worse. However, 
the complete Palomar-WIRC catalog covers a total area 
of 1.47 deg 2 out of which EGS is only a small fraction 
(0.20 deg 2 ). At j^vEGA^20 our resu lt s are also consis- 
tent w ith t hose in | Hartlev et alJ (l2008f).lMcCracken et al.l 
(|2010f) and lForster Schreiber et all (|2006f) . following the 
same trend as the latter up to Xvega=21. From this 
comparison, we conclude that the IRAC catalog limited 
to [3.6]=23.75 is a good proxy of a if-limited sample with 
at least ifyEGA<21mag. 

The top-right and the bottom panels of Figure [TT] 
shows the comparison of the number counts for DRG and 
BzK (restricted to z>1.4) galaxies with other values from 
the literature. There is good agreement within the typi- 
cal scatter (0.1-0.2 dex), generally associated with cosmic 

1 We used the following transformations when required 
AK V bga-Ab(UKIRT,CFHT,SOFI)=1.90,1. 85,1.87. 
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variance. Our counts reproduce the most representative 
features of the overall distribution, namely: the plateau 
in DRG and p-BzK around i^vEGA^20.5 and the steep 
slope in the counts of s-BzK. We n ote that our s-BzK 
cou nts are slightly above th ose from lBlanc et al.l (2008) 
and lMcCracken et al.1 (|2010t ) which count with very large 
surveyed areas (0.71 deg 2 and 2deg 2 , respectively). On 
the contrary, o ur res ults are in excelle nt agreement with 
Hartley et all (|2008l 0.63deg 2 ). In IMcCracken et al.l 
(2010) the a uthors argue that th eir disagreement with 
the counts of lHartlev et al.l (|2008|) is the result of an in- 
correct color correction in the filter system. However, for 
this w ork we used the exact same fil ters as in Dadd i~et al.l 
( 2004) obtaining similar results to lHartlev et al.l . Thus, 
the most plausible explanation is that there is an excess 
of galaxies at z*~1.5 in our region. 

For p-BzK our resu lts are lie betw e en th ose of 
lHartlev et all (|2008l ) and lMcCracken et al.l(j2010D . How- 
ever, the counts of p-BzK exhibit the largest scatter of 
the three populations. This is not surprising given that 
p-BzK target a m ore constrained popu l ation, prone to 
stronger clustering (IGrazian et al.ll2006bllKaiisawa et al.l 
l200lHartlev et al.1 12001 IMcCracken et al.l 12010ft and 
hence significantly affected by LSS. Table [5] sum- 
marizes the accumulated surface densities of DRG 
and BzK galaxies up to i\~vEGA=20 and JCvega=21. 
The values are roughly consistent with the results 
of the studies shown in Figure [11] and with other 
values from the literature (0.89 DRG/arcmin 2 in 
IQuadri et all 120071: 3.1 s-BzK/arc min 2 and 0.24 p- 
BzK/arcmin 2 in iReddv et al.l I2005t3.2 s-Bz K /arcmin 2 , 
0.65 p-BzK /arcmin 2 in iGrazian et al.l l2007f) . As men- 
tioned above, the excess of s-BzK by a factor of ~1.5 
could be caused by a source over-density in the area. 
Note that this excess does not necessarily affect s-BzK 
and p-BzK in the same manner due to the different clus- 
tering properties of each population (lHartlev et al.|[2~008l 



IMcCracken et all 2010. Ha rtley et al.H2010D In fact, our 
density of p-BzK is not among the lowest values. 

Figure [T2~l shows the z p hot distribution of DRG, s-BzK 
and p-BzK galaxies with ifvEGA <20 and ifyEGA <21, 
compared to some results from the literature. The dis- 
tributions are convolved with a 5z=0.1 kernel in order to 
account for the ;z p hot uncertainties. The redshift range 
spanned by the different galaxy populations is in good 
agreement with the usual distributions observed in other 
studies, i.e., z>2 for DRG and 1.4<z< 2.5 for BzKs 
(|Daddi et al.l 12004). DRG present a secondary redshift 
peak around z~l, that accounts for a significant frac- 
tion of the total population at bright (ifyEGA<20) mag- 
nitudes (as already pointed out by other authors, e.g., 
IQuadri et aT] 120071 IConselice et al.il2007ft . Nevertheless, 
our surveyed area (0.35 deg 2 ) is not large enough to make 
(bright) low-z DRG the dominant fraction, as in the 
0.70 deg 2 of the if-band Palomar survey where ^70% o f 
these galaxies are found at z<1.4 [C~onselice et al.l (|2007j ). 
As expected, s-BzK and p-BzK present almost identical 
redshift distributions, although the latter seems to have 
a more extended high-redshift tail, being also less prone 
to low redshift interlopers (probably as a consequence 
of the more restrictive color criteria). In summary, our 
results about the number density and redshift distribu- 
tion of color-selected z >l galaxy samples are consistent 
with previous studies (jReddv et al.ll2005l IGrazian et al.l 
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Fig. 12. — Photometric redshift distributions of s-BzK, p-BzK 
and DRG galaxies (from top to bottom) drawn from the IRAC- 
selected sample with [3. 6] < 23. 75. The black line and gray line show 
the distributions at _R"vEGA<21 and XyEGA<20, respectively. 
Our results are compare d with the median ( and q uartile) of the 
redshi ft distributi ons in | Gr azian et al. (20Q7], red), IQuadri et alj 
J2007I . green) and IMcCracken et alj J2010I . blue) (also top to bot- 
tom). The two (green) intervals in the redshift distribution of DRG 
indicate the median values of the distribution at redshifts lower and 
higher than z=1.5 in the work[Quadri et al. (2007]). 



TABLE 5 

Surface density of DRG and BzKs 





^VEGA<21 


^VEGA<20 






z h 


p a z° 


DRG 


1.4 


2.47 


0.5 2.24 


s-BzK 


5.0 


1.89 


1.5 1.70 


p-BzK 


0.5 


1.85 


0.3 1.73 



Note. — 

a Surface density of DRG and BzK in arcmin -2 . 
b Median photometric redshift of each sub-sample. 

l2007t PG08), indicating that the photometric redshift 
estimates are generally robust at high-redshift. 

3.7. Photometric redshift distribution 

Figure [TBI shows the z p hot distribution for the IRAC se- 
lected sample in the main region, limited to [3.6]<23.75. 
In addition, we also plot the redshift distributions of the 
galaxies detected at 24 /xm, 70 /im and the sub-sample 
with spectroscopic redshifts. In order to derive a re- 
alistic distribution, accounting for the uncertainties in 
Zphot, the distribution was convolved with the typical 
width of the zPDFs. We used a conservative upper limit 
of Az/(1 + z)=0.07. The shape of the distribution is 
consistent with that expected for a magnitude limited 
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Fig. 13. — Photometric redshift distribution for the IRAC 
([3.6]<23.75; black), MIPS 24 /im (x3, red), MIPS 70 fan (x6, 
orange) and spectroscopic (grey area) samples in the main region 
of the EGS. The distribution of z phot in 106 («' <24.5) for the 
overlapping area with the IRAC mosaic is shown in blue for com- 
parison. 

sample. At low redshift the number density increases as 
we probe larger volumes, and then an exponential decay 
is observed as the sources get fainter and the detection 
probability decreases. 

The position of the minor prominences in the z p hot dis- 
tribution are roughly consistent with the most remark- 
able peaks observed in the spectroscopic redshift dis- 
tribution at z^0.3, z^0.7, z~l. The median redshift 
of the photometric redshift distribution is z—1.2, 75% 
of the sources are below z=2.1, and 90% below z=2.7. 
The medi an of the d i stribu tion is consistent with the 
results of lllbert et all (120091) in the COSMOS field for 
an i'+3.6/im selected sample (i'<25, f(3.6)>l £tJy). Al- 
though the IRAC S-COSMOS catalog is shallower than 
ours, with a ^50% completeness level at f(3.6)=l /iJy, 
the median redshift limited to their faintest magnitude 
bin (24.5<z'<25), z =1.06, is similar to ours. Note that 
the small differences could arise from the presence of un- 
derlying LSS in EGS, whereas this effect is largely re- 
duced in the COSMOS sample due the larger area of 
field (^1.73 deg 2 ). Finally, the distribution is also in 
good agreement with our results in PG08 for the averaged 
redshift distribution of a combination of IRAC-selected 
catalogs in the HDFN, CDFS and Lockman Hole fields. 
The total combined area in PG08 is approximately that 
of the main region of EGS, and the limiting magnitude 
of the catalog was slightly lower (f(3.6)<1.6/xJy). How- 
ever, the redshift distribution also peaks around z=0.8-l, 
consistently with ours. 

4. SED ANALYSIS: STELLAR MASSES 

In this section we describe the method used to estimate 
stellar masses based on the SED fitting. In addition, 
we analyze the goodness of our stellar mass estimations 
quantifying the systematic and random errors linked to 
assumptions in the input parameters for the stellar pop- 
ulation modeling. For the discussion in § 14.2.11 we use 
only the spectroscopic sample in the main region (which 
count with better photometry) and we force the z p hot to 
the spectroscopic value. 



4.1. Stellar mass estimates 

The stellar mass of each galaxy is estimated from the 
wavelength-averaged scale factor required to match the 
template monochromatic luminosities to the observed 
fluxes. This is possible because our templates are ob- 
tained from stellar population synthesis models which are 
expressed in energy density per stellar mass unit. Note 
that the stellar mass estimate is not obtained from a 
single rest-frame luminosity and its corresponding mass- 
to-light ratio, which has been a typical procedure seen in 
the literature, but from the whole SED. In our method, 
the fit to the multi-band data implicitly constraints the 
mass-to-light ratio by determining the most suitable tem- 
plate. Then, we estimate the mass from the averaged 
template normalization, weighted with the photometric 
errors. This approach is less sensitive to the effects of 
the star-formation history (SFH) or the photometric and 
template uncertainties in a single band. Objects fitted 
with pure AGN templates have no stellar mass estimate, 
as their SED is dominated by non-stellar emission. The 
random uncertainty of the stellar mass is estimated with 
a bootstrap method by randomly varying the photomet- 
ric redshift and observed fluxes based on their quoted 
errors. 

4.2. Accuracy of the stellar masses 

In addition to the uncertainties inherited from the 
probabilistic nature of z p hot and the intrinsic photomet- 
ric errors, there is another source of systematic uncer- 
tainty associated to the assumptions in the SED model- 
ing. Although significant effort has gone into providing 
accurate SPS models, key ingredients of the theoretical 
predictions are still poorly understood. As a result, there 
can be substantial differences in the physical properties 
estimated with many of the well-tested SPS models avail- 
able in the literature. Most of these differences arise from 
the different parametrizations of potentially uncertain 
phases of the stellar evolution, such us the asymptotic 
giant branch (AGB) or the thermally pulsating AGB 
(iMarastonl 120051 iBruzuall I2007L iKannappan fc Gawiserl 
120071) Another critical aspect, is the choice of an IMF. Al- 
though this is essentially assumed to introduce a change 
in the overall normalization of the stellar mass, there are 
additional effects attached, e.g., a change in the balance 
between low-mass and high-mass stars varies the rela- 
tive fraction of stars in different points of the isochrones. 
Thus modifying the colors and M/L of the modeled 
galaxies at different evol utionary stages (|Marastonlll998l 
Ivan Dokkum et al.ll2008l) . 

Apart from the choice of SPS models and IMF, ad- 
ditional effects might arise from the assumed SFH, 
usually parametrized with r-models, or the choice of 
a dust e xtinction law and metall icity. As recen tly 
shown in iMuzzin et all ()2009f ) (also iKriek et alJl2008l or 
iMarchesini et al.ll2009|) . due to all these effects, the phys- 
ical properties of galaxies estimated from broad-band 
photometry often presents large uncertainties (typically 
within a scatter of 0.2 dex for stellar masses), in addition 
to systematic offsets. Moreover, these uncertainties can 
be even larger (up to 0.6 dex) for particularly sensitive 
galaxy populations at high-z, such as bri ght red galax- 
ies. See for e xample the series of pap er by IConrov et al.l 
(f2009l l2010h : IConrov fc Gum] §010) for a detailed dis- 
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Fig. 14. — Left: Comparison of the stellar masses estimated with [P01,SALP,CAL01] using spectroscopic and photometric redshifts in the 
main (black) and flanking (red) regions, respectively. The blue line indicates the best Gaussian fit to the central values of the main region 
distribution. Right: Stellar masses of the galaxies in the sample as a function of redshift. The black dots depict galaxies with [3.6]<23.75 
(85% completeness level of the sample). The red dots depict galaxies 23.75<[3.6]<24.75 (3cr limiting magnitude). The green dots show 
galaxies with spectroscopic redshifts. The blue lines indicate the 90% and 10% percentiles of the mass distribution as a function of redshift 
for the galaxies with [3.6]<23.75 (solid) and 23.75<[3.6]<24.75 (dashed), respectively. 



TABLE 6 

Comparison of stellar masses computed with different 
modeling assumptions 



IMFs 

(1) 


SPS model 

(2) 


Dust 
(3) 


A log(M) 
(4) 


Alog(M) 
(5) 


SALP-KROU 
SALP-KROU 
KROU-CHAB 


CB09 
P01 
CB09 


CAL01 
CAL01 
CAL01 


0.19±0;0V 
0.03± ; 
0-04±«; f ± i 


0.19±0;1| 
0.13±S-I? 

o.o7±»;?i 


SPS model 


IMF 


Dust 


A log(M) 


Alog(M) 


BC03-CB09 
P01-CB09 
P01-M05 
CB09-M05 


CHAB 
KROU 
KROU 
KROU 


CAL01 
CAL01 
CAL01 
CAL01 


0.04±°-ff 
0.15±° : | 
0.39±R : || 


0.07±S|? 
0.08±g'i 
0.30±S : p 
0-20±°l 


Dust 


IMF 


SPS model 


A log(M) 


Alog(M) 


CAL01-CF00 


SALP 


P01 


-0.03±°;i 


0.00±°;§° 



Note. — Comparison of the stellar masses obtained under different 
combinations of the modeling assumptions. 

(1),(2),(3) SPS model, IMF and dust extinction law, alternatively. 
The first column indicate the parameters being compared. 

(4) Log of median value and quartiles of the difference for galaxies 
with log(M)<lOAl0. 

(5) Same as (4) for galaxies with log(M)>10A4o 



cussion of all these issues. 

Taking these considerations into account, in the fol- 
lowing Sections, we analyze the accuracy of our stellar 
mass estimates quantifying the uncertainty budgets as- 
sociated with different effects. First, we study the effect 
of photometric redshift uncertainties. Then, we evalu- 
ate the impact from the choice of SPS models, IMF and 
dust extinction law restricting the analysis to the spec- 
troscopic sample. For the sake of clarity, we refer all 
comparisons to a default choice of SED modeling param- 



eters (as described in [3]) characterized by SPS models, 
IMF and extinction law [P01,SALP,CAL01]. Finally, we 
verify that our stellar masses provide realistic values by 
comparing them to other stellar mass catalogs available 
in the literature. 

Note that although the a priori assumptions on the 
SFH can also introduce systematic effects in the esti- 
mated stellar masses, an in depth analysis of these is- 
sues is clearly beyon d the scope of this paper (see e.g., 
iMaraston et alJfeOlOl for a detailed discussion) . Nonethe- 
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Fig. 15. — Comparison of the stellar masses obtained using different IMFS and dust extinction laws for a given SPS library. Top-left: 
P01 models with a SALP and KROU IMFs. Top-right: CB09 models with a SALP and KROU IMFs. Bottom: CB09 models with a KROU 
and CHAB IMFs. Bottom-right: P01 models with a CAL01 and CF00 extinction laws. The histograms in the right part of the plot depict 
the ratio of the stellar masses obtained with each IMF for galaxies with M mot j e i<10A4Q (empty) and M mot j e i>10A'fQ (filled). The solid 
lines above the histogram show the median value and la of the distribution at both sides of the median. 

less, a comparison of the results obtained with a single 
exponentially-declining stellar population (1-POP) and 
with a single population plus a second burst (2-POP) is 
presented in PG08 (Appendix B) along with similar tests 
to the ones presented in the next section. 

The catalog of stellar masses presented in this paper 
(see § O contains the different values obtained with all 
the modeling configurations discussed in the next sec- 
tions. 



4.2.1. Effects of the Photometric redshifts, SPS models, 
IMF and extinction law 

The left panel of Figure[T4lshows the scatter in the stel- 
lar masses estimated using z p hot and z sp0 c for the 7,636 
spectroscopic galaxies in the main (black) and flank- 
ing regions (red). Approximately 68% and 90% of the 
sources are confined within a rms of 0.16, 0.34 dex and 
0.20, 0.39 dex in each region, respectively. Nonetheless, 
the distribution shows a pronounced central peak that 
it is well reproduced by a Gaussian distribution (blue 
line) with extended wings, indicating that for the most 
accurate redshifts, the scatter is substantially reduced 
(^0.065 dex). This is in good agreement with the re- 



sults of lllbert et all (|2010l see Figure 3) scaled to the 
overall accuracy of our photometric redshifts, which is 
slightly lower. The right panel of Figure [14] shows the 
range of stellar masses as a function of rcdshift for the 
whole sample (black). In order to illustrate the ap- 
proximate limiting stellar mass inherited from the mag- 
nitude limit ([3.6]<23.75; 85% completeness), we also 
depict the galaxies up to the 3cr limiting magnitude 
(23.75<[3.6]<24.75, red dots). Approximately 90% of 
the galaxies with [3.6]<23.75 present log(M)>10M s at 
z>2.5 (blue line) in agreement with our results in PG08 
for a similar limiting magnitude. Similarly, ~10% of the 
faintest galaxies (23.75<[3.6]<24.75), absent in our main 
sample, present stellar masses larger than log(M)>10Af Q 
(blue dashed line). Note however that the completeness 
in stellar mass can not be directly extraploted from these 
limits because for any given redshift, galaxies with dif- 
ferent ages present different mass-to-light ratios. Hence, 
the completeness is an age (or color) dependent value. 
In particular, magnitude limited samples are known to 
be incomplete again st the oldest (red) galaxies (see e.g., 
iFontana et al.l l2006f) . A detailed analysis of the com- 
pleteness limit as a function of the galaxy type will be 
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Fig. 16. — Comparison of the stellar masses obtained using different SPS models with the same IMF. Top-left: BC03 and CB09 models 
with a CHAB IMF. Top-right: P01 and CB09 models with a KROU IMF. Bottom-left: P01 and M05 models with a KROU IMF. Bottom- 
right: CB09 and M05 models with a KROU IMF. The histograms in the right part of the plot depict the ratio of the stellar masses obtained 
with each model for galaxies with M mo ^ B i<10M<3 (empty) and M mo( j o i>10.Mo (filled). The solid lines above the histogram show the 
median value and lcr of the distribution at both sides of the median. 



included in a forthcoming paper. 

The first test on the effect of the SED modeling as- 
sumptions consist on a comparison of the stellar masses 
compute d with t hree differe nt choice s of th e IMF; a 
SALP , IKroupal (|200ll) and IChabrierl (120031 hereafter 
KROU,CHAB) IMFs. The naive expectation is that, the 
stellar masses obtained with a SALP IMF are on aver- 
age larger than those obtained with the other two, as it 
predicts a larger number of low- mass stars. On the con- 
trary, the IMFs of KROU and CHAB are quantitatively 
very similar and therefore the differences are expected 
to be small and mass independent. Table [5] summarizes 
the median value and quartiles of the comparison of stel- 
lar masses obtained with each IMF in combination with 
the P01 and Chariot & Bruzual (2009; CB09) models 
and a CAL01 extinction law against the reference val- 
ues. The top panels in Figure fTS] show this comparison 
for SALP or KROU IMFs (left) and a KROU or CHAB 
IMFs (right) and the CB09 models. In both cases the dif- 
ference is essentially a constant value of factor ~1.6 and 
1.2, respectiv ely, consistently with the results of the lit- 
erature (e.g., [S alimbcni et al.||2009|. iMuzzin et aT1l2009L 
iMarchesini et al.l 12009). On the contrary, the difference 



in the values obtained with a SALP or KROU IMFs for 
the P01 models is significantly smaller than the for the 
CB09 models, showing also a larger scatter and a depen- 
dence on the stellar mass. This effect seems to be related 
with a difference in the age dependency of the mass-to- 
light ratio_Jor ea£hIMF__^ (see 
e.g., lMaraston|[l998l Ivan Dokkum et al.l [20081 for a de- 
scription of these effects). 

The second test on the modeling parameters is the 
comp arison of stellar ma s ses ob tained with the models of 
POl. lBruzual fc Charlotl (12001 . BCO.Ti. iMarastoi] ([2005L 
M05) and CB09 for the same IMFs. The models of M05 
were the first to account for the contribution of the TP- 
AGB phase in the SPSs, a consideration that is expected 
to lead to lower stellar masses compared to those ob- 
tained with P01 and BC03. However, this difference 
should be reduced in the CB09 models, the updated ver- 
sion of BC03, which include an improved treatment of 
this particular phase. The overall results of the compar- 
ison between models is also summarized in Table [5] and 
in the panels of Figure [TBI 

Interestingly, we find that the difference between 
[BC03,CHAB] and [CB09,CHAB] is relatively small, 
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Fig. 17. — Difference between our bes t-fit stellar mass usin g [P01 
mass limited sample (log M>11X ) of ITruiillo et al.l H2007T) . The 
the comparison to Trujillo et al. is limited to log M>11 Mq. 

~0.04dex, and mostly independent of the stellar mass. 
This suggests that, at least for the present sample, tak- 
ing into account the TP-AGB phase does not introduce 
significant differences. A possible explanation could be 
that, since the spectroscopic sample consist mostly on 
z<l galaxies, the available photometric coverage is not 
probing the rest-frame NIR with sufficient detail. Only 
at higher redshifts (z>l) the IRAC bands would start 
probing the region of the SED that heavily affected by 
the TP-AGB phases. Note also that the spectroscopic 
sample analyzed here might not be a critical population 
to constraint the effect of the TP-AGB, as for exam- 
ple th e post-starburst galaxies studied in IConrov et al.l 
(|2OT0l) . 

In addition, we find that the estimates with 
[P01,KROU] are larger than those obtained with 
[CB09,KROU] and [M05,KROU] with an average offset 
of 0.15 dex and 0.39 dex, respectively. The difference with 
respect to M05 is consistent with previous res u lts (e.g., 
Maraston etsll 120061 Ivan der Wei et al.l 120061 iBruzuall 
2007|) in spite of the slight dependence on the mass. 
However, the 0.16 dex offset between [CB09,KROU] and 
[M05,KROU] (illustrated for completeness in bottom- 
right panel of Figure [16]) is larger than expected reveal- 
ing a more complex relative difference between the two 
libraries beyond the treatment of the TP-AGB phase. 

Finally, the bottom-right panel of Figure [15] shows the 
comparison of the stellar masses estimates obtained with 
a CAL01 and a Charlot&Fall (2000; CF00) dust extinc- 
tion laws for the P01 models and a KROU IMF. The 
most relevant differences between a CAL01 and CF00 
extinction laws is that the latter presents a larger atten- 
uation of the stellar component, which effectively leads to 
lower fluxes (mostly in the UV) for similar values of the 
extinction. Furthermore, the wavelength dependence of 
the attenuation in CF00 is greyer (i.e., shallower) than 
in CAL01. The overall result of the comparison is an 
small offset of -0.03 dex with a ~0.2 dex rms, similar to 
what we found in PG08. This is also in good agreement 
with the results of iMuzzin et al.l (j2009f ). indicating that 
the treatment of the extinction law does not play a ma- 
jor role in the estimate of the stellar mass (although it is 
more relevant for other estimated parameters). 

In summary, we find that after accounting for the dif- 
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,SALP] and the stellar masses of lBundv et al.l j200Cf;). (lef t), and the 
histograms in the right side are the same as in Figure [151 Note that 

ferent systematic offsets, all models seem to be roughly 
consistent within a factor 2 (~0.3dex). However, there 
are mass dependent systematics that should be taken into 
account in th e analysis of overall properties of galaxy 
samples (e.g.. iMarchesini et af1l2009L iTavlor et alJl2009l 
lllbert et al.lboiOD . 

4.2.2. Comparison to other stellar mass catalogs 

In this section we compare our st ellar masses with 
the es timates from lBundv et al.l (]2006t ) and ITruiillo et al.l 
(2007). In the former, the authors derived stellar masses 
for a large sample of galaxies with spectroscopic red- 
shifts from the DEEP2 survey in the EGS. In the lat- 
ter, the authors combined spectroscopic and photometric 
redshifts to study the properties of a mass limited sample 
(log M>11 yVf©). Both works used the same photomet- 
ric dataset consisting on 5 bands: BRI from the CFHT 
survey, and JK from the Palomar NIR survey. The stel- 
lar masses in both cases were essentially computed based 
on the fitting of the SEDs to a grid of tem plates derived 
from BC03 models with a lChabrierl (]2003t) IMF and ex- 
ponentially decreasing SFHs. In particular, Bundy et al. 
used the rest-frame if-band luminosity and mass-to-light 
ratio to scale the templates and compute the probability 
distribution of the stellar mass and the most likely value. 
On the contrary, Trujillo et al. (based on the results of 
Conselice et al. 2007) did not renormalize the templates 
in a single band but used the whole SED to scale the 
fluxes, similarly to our approach but restricted to only 
5-bands. 

We cross-correlate the catalogs using a 2" radius, and 
we double check the validity of the match ensuring that 
the spectroscopic redshifts (independently matched) are 
the same. The final sample contains 4706 and 791 galax- 
ies detected in the catalogs of Bundy et al. and Tru- 
jillo et al., respectively. For the comparison to Bundy et 
al. the photometric redshifts were forced to the spectro- 
scopic value and for the comparison to Trujillo et al. the 
photometric redshifts were forced to the values quoted in 
their paper. 

The left panel of Figure [T7] shows the comparison of the 
stellar masses with Bundy at al. for our default model- 
ing assumptions. Our estimates are slightly lower with an 
median difference of A log(M)=-0.07±0.21 dex. Also, we 
find that the stellar masses computed with [P01,SALP] 



are in better agreement with Bundy et al. than those ob- 
tained using the same modeling configuration as in their 
work, [BC03, CHAB], which would increase the differ- 
ence in smaller masses to Alog(M)= -0.12 dex. We fur- 
ther investigate if this offset is caused by a difference in 
the photometry by comparing our if-band magnitudes to 
those of Bundy et al. that were computed using 2"radius 
apertures (for the SED fitting) . The sources in Bundy et 
al. are on average AK =0.12 mag fainter than in our cat- 
alog, which would imply a larger difference in the stellar 
masses if we simply scale their magnitudes to our pho- 
tometry. Thus, the most plausible explanation for this 
small offset is the use of different techniques for estimat- 
ing the stellar masses, and specifically the use of IRAC 
data in our study. The right panel of Figure [IT] shows 
the comparison of the stellar masses with Trujillo at al.. 
The overall comparison presents a good agreement with 
a median difference of Alog(M)=0.10±0.25dex, slightly 
larger than the offset to masses of Bundy et al. However, 
the scatter of the distribution is quite similar to that of 
the comparison to Bundy et al. for the highest stellar 
masses log(M)>ll M Q . 

5. SED ANALYSIS: STAR FORMATION RATES 

In this section, we present the estimations of the SFRs 
of the galaxies in our IRAC sample based on their UV-to- 
FIR SEDs. We also discuss the quality of these estimates 
as well as their associated systematic uncertainties. 

The SFR of a galaxy is frequently computed from the 
UV and IR luminosities through theoretical or empiri- 
cal calibrations. As young stellar populations emit pre- 
dominantly in the UV, this wavelength range is highly 
sensitive to recent events of star formation. However, 
this UV emission is usually attenuated by dust, which 
re-emits the absorbed energy in the thermal IR. Con- 
sequently, the ongoing SFR can be estimated either by 
correcting the UV luminosity for extinction or combining 
the IR emission and the unobscured UV flux. 

Here we focus on the latter approach making use of the 
high quality FIR fluxes observed with MIPS at 24 and 
70 /urn. Thus, assuming that the total SFR of a galaxy 
can be estimat ed by summing up two comp onents (see 
e.g.. iBell et "all 120051; Uglesias-Paramo et al.ll2007t ): the 
part of the star formation that is probed by a tracer af- 
fected by dust attenuation, so we only are able to observe 
directly a fraction of it (i.e., the unobscured component), 
and the part of the star formation that is hidden by dust 
(obscured component). The unobscured star formation 
can be measured with the rest-frame UV emission, which 
can be estimated from the optical/NIR SEDs for the 
galaxies in our sample. The obscured component can 
be estimated from the total IR thermal emission (thus, 
we will refer to it as IR-SFR or IR-based SFR). How- 
ever, its calculation is usually affected by the choice of 
template libraries fitting the IR part of the SED and, 
more significantly, by the photometric coverage in the 
MIR-to-mm spectral range. In this section, we will fo- 
cus on the analysis of the IR-based SFR and the random 
and systematic uncertainties associated with the differ- 
ent procedures used to estimate it. 

The structure of this section is as follows. First, we 
describe how we fit the IR part of the SEDs to dust emis- 
sion models, and present the different methods used to 
estimate and IR-based SFR from monochromatic and in- 
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Fig. 18. — Ratio of the two components of the total SFR 
(SFR=SFRj7v r ,o6s+SFR/^j; Equation |4j as a function of the 
SFRtir for a sub-sample galaxies detected in MIPS 24 fim 
(f(24)>60/iJy).Each contour contains (from the inside out) 25%, 
50%, 75% and 90% of the sample, respectively. The underlying 
black dots depict the individual values of the ratio of SFRs. The 
grey stars with error bars depict the median and lcr of the ratio of 
SFRs in bins of SFRtir. 

tegratcd luminosities in the MIR-to-mm range. Then, we 
compare these different methods and discuss the system- 
atic and random uncertainties inherent to the calculation 
of IR-based SFRs. 

In this section, the reader must have in mind that 
the most useful information to estimate IR-based SFRs 
comes from the MIPS 24 /im fluxes. The reason is sim- 
ple: these observations are the deepest in the MIR-to- 
mm range, so we only have this SFR tracer for the vast 
majority of sources in our sample. Ideally, it would be 
desirable to have other fluxes in the IR to constrain the 
fits to dust emission models, but this is only possible for 
a very small fraction of galaxies which have MIPS 70 /im 
data, or other photometric points in the (sub)-mm. Even 
with Herschel data, there will be a significant population 
of galaxies that will only count with the MIPS 24 /im 
flux. Therefore, an important part of our discussion will 
be assessing the reliability of IR-based SFRs based only 
on MIPS 24 fj,m data. For that purpose, we will take 
advantage of the very deep observations carried out at 
70 fxm within the Spitzer FIDEL Legacy Project, study- 
ing the variations in the estimated IR-based SFRs fitting 
MIPS 24 /mi and MIPS 70 /jm simultaneously. 

5.1. IR SED fitting 

Typically, IR-based SFR are computed either from the 
total IR luminosity, i.e., the integrated emission from 8 
to 1000 /im [L(TIR)], or from monochromatic luminosi- 
ties at different wavelengths. Both methods require a 
detailed characterization of the IR SED, which is usu- 
ally obtained by fitting the observed fluxes to dust emis- 
sion templates. However, as mentioned above, these es- 
timates are largely dependent on the choice of templates. 
An issue that is usually aggravated by the fact that typ- 
ically the only measurement of the MIR emission comes 
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from the 24 /an data, and occasionally 70 /im, whereas 
the total IR luminosity is commonly dominated by the 
emission at longer wavelengths A^100/im. 

Thus, in order to study in detail the intrinsic uncer- 
tainties in the IR-based SFRs arising from these issues, 
we follow two different approaches to fit the IR data 
to the dust templates: I) we study the galaxies de- 
te cted at MIPS 24 fitting only this flux to models 
oflCharv fc Elbazl (l200ll CEOf hereafter). iDale fc Heloul 
d2002j, DH02 hereafter), which is a usual scenario in stud- 
iesof the IR-emiss i on at high-z (see e.g. ISantini et al.l 
120091 IWuvts et all 120081 ). In this case, we asses the 
differences between IR-based SFR (hereafter SFR, (24)) 
estimated with several methods, and the impact of 
using different models; 2) we restrict the analysis to 
galaxies simultaneously detected in IRAC and MIPS 
24 and 70 /im, fitting all fluxes at rest-frame wave- 
lengths A> 5 /im (where the luminosity of a galaxy must 
present a significant non-stellar contribution; see, e.g., 
IPerez-Gonzalez et all [20061 ISmith et all 120071 1 to the 
models of CE01, DH02 and also lRieke et al.l (|2009l R09 
hereafter). We refer to these galaxies as the best-effort 
sample and their SFRs(8,24,70). The notation indicates 
that the fit essentially includes 8, 24 and 70 //m data 
up to z^0.6, and 24 and 70 /im data at higher redshift. 
Based on this sample we can study the impact of having a 
better constrained IR SEP against the M IPS 24 /im only 
scenario (e.g, as in lKartaltepe et aU ^lO). In both cases, 
the fitting is carried out by fixing the redshift to z p hot 
or z spcc (if available). Then, the excess resulting from 
subtracting the predicted contribution from the stellar 
flux (given by the best fitting optical template) to the 
MIR bands is fitted to each set of models. In the case 
when only MIPS 24 /im data is used, the templates are 
not fitted but rather scaled i.e., we obtain the rest- frame 
monochromatic luminosity for that flux and redshift and 
we select the most likely template based in their abso- 
lute normalization in the total IR lumin osity (as in e.g., 
IPapovich etall 120061 ISantini et al.ll2009l) . Moreover, for 
sources undetected in the 24 /im data, we set an upper 
limit of f(24)=60/i Jy, the approximate SNR=5 level of 
the MIPS data in EGS (see paper I for more details), 
which allows us to provide an upper limit of the IR-based 
SFR. Figure [2] shows and example of the IR SED fitting 
jointly with the optical template. 

Based on the best fitting templates, we computed sev- 
eral IR monochromatic and the integrated luminosities 
[L(A) and L(TIR), respectively] as the median value of 
all the fitted template sets. In the following, we describe 
various possibilities for IR-based SFRs based on L(A) at 
different wavelengths. This relations are calibrated from 
galaxy samples counting with extensive IR coverage (at 
least more than 3-4 bands), and provides an alternative 
estimate of the SFR based on milder template extrapo- 
lation than L(TIR), which in principle makes them more 
robust when only few bands are available for the fitting. 
Note however that the rest-frame wavelengths around 10- 
30 /jm are wildly variable and thus extrapolating lumi- 
nosities in this region involves significant uncertainties, 
e.g., L(8) based on 24 /im data at z^2 (see § 15.3.21 for 
more details). 



5.2. Total SFR and IR-based SFRs estimates 



Our method to estimate the total SFRs is based 
on a combination of the I R emission and the unob- 
scured UV flux (similarly to IPerez-Gonzalez et al.ll2006l 
iKennicutt et al.ll2009l . ISantini et al. |200£ ).In p articular, 
we use the prescrip tion of iBell et al.l ( 2005 . see also 



IPapovich et al.ll2007l ). which is based on th e calib ration 
for the total IR luminosity of IKennicutt JT998) , and 
parametrizes the contribution of radiation that escapes 
directly in the UV. 

SFR = SFRtir + SFR(jv,obs (3) 

SF^Moyr- 1 ) = 1.8x I0" 10 [L(T/i?)+3.3x L(0.28)]/i Q 

(4) 

where L(TIR) is the integrated total IR luminosity, and 
L(0.28) is the rest-frame monochromatic luminosity at 
0.28 /im (uncorrected for extinction). The well sampled 
SEDs of our galaxies at optical wavelengths allow a ro- 
bust estimation of L(0.28) by interpolating in the best 
fitting optical template. However, as described in the 
previous section, the value of L(TIR) is strongly model 
dependent, as it is based on an extrapolation from one or 
a few MIR fluxes to the total emission from 8 to 1000 /im. 
An alternate possibility is to obtain other IR-based SFRs 
based on L(A), thus reducing the template dependence. 
In the following, we will refer to the IR-based SFR de- 
rived from L(TIR) as SFRtir. In addition, we compute 
four other IR-based estimates. 

The first estimate is based on rest-frame monochro- 
matic luminosity at 8 /im (hereafter SFRbos)- These esti- 
mate make use of the em pirical relation between L(8/im) 
and L(TIR) described in lBavouzet et al.l (|2008l ) and the 
Kennicutt factor to transform to SFR. 

SFR BQ8 (M Q yr- 1 ) = 1.8 x 10~ 10 x (377.9 x i(8)°- 83 )/i Q 

(5) 

The seco n d met hod in based on equation (14) of 
iRieke et al.l (|2009l ). that relates the SFR (hereafter 
SFR RO g) to the observed flux in the MIPS 24 /xm band 
and the redshift. The redshift dependent coefficients of 
the relation were computed using averaged templates de- 
rived from a set of empirical IR-SEDs fitting local galax- 
ies. This estimation of the SFR is independent of the 
rest, as it is based on different templates. The conversion 
from IR luminosities to SFRs is also computed through 
the Kennicutt factor. However, the authors scaled the 
factor to a Kroupa-like (2002) IMF (the original factor is 
for a SALP IMF) mu ltiplying it by 0.66 (a similar con- 
version is obtained in lSalim et al.ll2007l ). Here we undo 
that change for consistency with the other methods that 
are computed using the factor for a SALP IMF. 

The last method is not strictly an IR-based SFRs but 
an estimate the gl obal SFR. It is based on th e empiri- 
cal relation given in lAlonso-Herrero et al.l (|2006l ) between 
the rest-frame monochromatic luminosity at 24 /im and 
the SFR (hereafter, SFRa-hog), 



SFR A _ m6 (M Q yr- 1 ) = 1.51 x 10~ 8 x L(24) 



0.871 



(6) 



This formula is based on the calibration of L(Paa) ver- 
sus L(24 /zm) obtained for a set of local ULIRGS using 
the Kennicutt (1998) relation bet ween L(Paa) an d SFR . 
A similar result was obtained by ICalzetti et all (|2007l ) 



25 



for resolved star-forming regions in local starburst (see 
also lKennicutt et aLll2009h . This estimation refers to the 
global SFR, n ot the IR-based SFR, becau se the empiri- 
cal relation in lAlonso-Herrero et al.l ()2006l ) already takes 
into account the unobscured star formation (measured 
through the observed Paa emission) and the extinction 
correction [applied to calculate L(Paa) in that paper]. 

5.3. Accuracy of the IR-based SFRs 

In the following Sections we analyze the systematic 
uncertainties in the IR-based SFRs associated with the 
use different models and indicators, and also the number 
of photometric bands available for IR SED fitting. 

First, we compare the values obtained with each of 
the methods presented in the previous section for a 
sub-sample of MIPS 24 /mi detected galaxies. In this 
case, the SED is fitted to 24 /tm data only (§ I5.3.1[) . 
Note that we have chosen several methods for estimat- 
ing the IR-based SFR that present intrinsically different 
approaches, using either integrated and monochromatic 
luminosities or observed fluxes. Here we also test the 
differences introduced by the use of the CE01, DH02 or 
R09 models. For simplicity, in this case, the comparison 
to the R09 models is done through the SFRs obtained 
with their empirical relation (SFRR09) instead of fitting 
the data to the three models. 

Second, we study the differences in the SFRs obtained 
for MIPS 24 /im sample and the best-effort sample, which 
count with a better IR coverage based on IRAC-8.0 plus 
MIPS 24 and 70 /im data (§ 15X2"]) . With this test we 
quantify the systematic effects associated with use of lim- 
ited IR data. Finally, we repeat the comparison of values 
obtained with each method for the best-effort sample in- 
cluding also highly accurate SFRs drawn from other au- 
thors based on a more detailed IR coverage (§ I5.3.3[) . 
Based on this comparison we asses the goodness of our 
best-effort SFRs and the reliability of the different meth- 
ods studied here. 

For the sake of clarity, we will refer all the compar- 
isons between the SFRs estimated with each method to 
SFRtir which, as explained in § 15 - H is computed from 
the average total infrared luminosity of all the fitted tem- 
plate sets. In addition, we will refer to them just as 
SFRs (dropping the IR prefix). In the case of SFRa-ho6- 
the proper IR-based SFRs is obtained by subtracting 
the contribution of the SFR[/y j0 b s - Nevertheless, as our 
working samples are composed by strong IR-emitters, we 
are biased towards dust obscured galaxies where this con- 
tribution is presumably small. For example, Figure 1181 
which shows the ratio SFHuv.obs / SFRtir, indicates that 
SFR[/y j0 fc s is lower than SFRtir (in most cases), with a 
clear trend for galaxies with intense star formation to 
present more and more extincted starbursts. 

5.3.1. Analysts of IR-based SFRs: MIPS 24 A* m sample 

Figure shows the comparison of the IR-based SFRs 
obtained with the different methods presented in Sec- 
tion 15.21 with respect to SFRtir as a function of SFR 
and redshift. All the estimates discussed in this sec- 
tion are based on 24 /im data only, i.e., SFRi(24). We 
omit the parentheses for simplicity. The SFRtir esti- 
mated separately with the models of CE01 (green dots) 
and DH02 (blue dots) are shown jointly with the me- 
dian and rms of both values in several SFR and lumi- 



nosity bins. The typical scatter of the SFRs estimated 
with both libraries is smaller tha n ~0.3 dex, consistent 
with the results by o ther authors (|Marcillac et al.l [20061 : 
iPapovich et~aLll2007D . When we compare the CE01 and 
DH02 libraries as a function of redshift, the maximum 
differences are observed for galaxies at z=l-2. For these 
sources, the estimates with the DH02 mo dels are larger 
than t hose with CE01 models, as found bv lSantini et al.l 
(2009). In this redshift range, the 24 /im band is probing 
the spectral region where the 9.6 /tm silicate absorption is 
found, jointly with the prominent PAHs around 7-9 /im. 
The shape of the models in CE01 and DH02 template sets 
is very different in this region, with the former presenting 
less prominent PAH features than the latter. In fact, all 
the DH02 models are identical below ^9 /im, while CE01 
models present a wide variety of spectral shapes, with a 
rising warm-dust continuum hiding the PAH features be- 
tween 6 and 20 /tm as we move to models with higher IR 
luminosities. 

The values of SFRa-ho6 are systematically smaller 
than SFRtir. For SFRtir > 100 M Q yr -1 , we find 
ASFR=-0.18±0.05dex. For smaller values of the SFR, 
where the unobscured and obscured star formation are 
comparable, SFRa-ho6 is down to a factor of 0.6 dex 
smaller than SFRtir, with a larger scatter. 

The comparison of SFRtir and SFRros clearly in- 
dicates that the em pirical relations L[8]-to-LTiR in 
iBavouzet et al.l (|2008l ) (Equation [SJ a nd in the models o f 
CE01 and DH02 (see e.g., Figure 8 of iDaddi et al|[2007]) 
are substantially different. The ratio of the two SFRs as 
a function of SFRtir is tilted with respect to the unity 
line, and consequently, both estimates are only consis- 
tent within a narrow interval around SFR^20 M Q yr _1 
(or z~l). For SFR T ir>100 and 1000 Moyr^ 1 (the lat- 
ter being the typical value for the z~2 galaxies detected 
by MIPS) the SFR B0 8 values are 0.3 and 0.8 dex lower 
than the SFRtir estimates, respectively. In contrast, for 
SFRTiR^lOMoyr" 1 , SFR B08 is larger than SFRtir by 
>0.2 dex. 

The equation to calculate SFRro9 (Ri eke et al.l 12009) 
varies with redshift. Consequently, the SFRro9 /SFRtir 
ratio presents different trends as a function of both lu- 
minosity and redshift. In terms of redshift, we distin- 
guish three regions: 0<z<1.4, 1.4<z<1.75, and z>1.75. 
At z=0-1.4, the ratio increases with redshift from an 
average value of -0.5 dex at z=0 to 0.5 dex at z^l.4, 
being close to unity at z^0.75. In the interval from 
1.4<z<1.75, the ratio decreases from 0.5 dex to nearly 
^0. Finally, at z>1.75, SFRr 9 values become roughly 
consistent with SFRtir with little scatter up to z=3, 
ASFR=0.02qq| dex. These large differences are related 
to the distinct shapes of the R09 and CE01/DH02 tem- 
plates. At z<0.5, the 24 /tm band probes a spectral range 
dominated by warm dust and emission features found by 
Spitze r at A~17 /tm and id entified with PAH or nanopar- 
ticles (| Werner et al.|[2004T ). At these redshifts, our sam- 
ple is dominated by galaxies with L(TIR)^10 10 L Q , and 
the CE01 models for these luminosities differ from the 
corresponding R09 templates by up to 0.5 dex in the 
A=16-20/im. This explains the differences at low red- 
shift in the right panel of Figure fT9l At z=0. 5-1.0, our 
sample is dominated by LIRGs, and CE01 and R09 mod- 
els for L(TIR)-1O 1O 75 L and L(TIR)-10 n ' 25 L are 
very similar (up to A=1.5mm), resulting on very similar 
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Fig. 19. — Comparison of different IR-based SFR indicators with respect to SFRtir as a function of SFRtir (left) and redshift (right) for 
galaxies detected in MIPS 24 (im (f(24 fim)>60fj,Jy). SFRtir is computed from the average value of Ltir in the t emplates of CE01 and 
DH02 fitted to the flux at 24 /im. The magenta points show SFRa— H06 estimated from the L(24) using the relation of Alonso-Hcrrero et al. 
(2006); the magenta stars and error bars in dicate the median value and la per SFR bin. The grey points show SFRrqs estimat ed from L(8) 
using the relation of Bavouzct ct al. (2008). The red points show SFRro9 estimated from MIPS 24 fj,m using the calibration of lRieke et alj 
(2009). The green and blue points depict SFRtir estimated from the templates of CE01 and DH02, respectively. The grey error bars 
depict the la uncertainty in SFRtir per SFR, bin. 



estimates of the SFR. At z~1.4, the 24 /im band probes 
the spectral region around 10 /im, and the galaxies de- 
tected by MIPS in this range have L(TIR)>10 115 L . 
For this luminosity, the CE01 and R09 models differ con- 
siderably due to the relative strength of the silicate ab- 
sorption. For example, for a L(TIR)>10 12 25 L Q , the R09 
template predicts a luminosity at 10 /im which is a factor 
of ^0.7 dex smaller than the CE01 model correspond- 
ing to the same L(TIR). Below 8 /im, rest- frame, the 
CE01 and R09 models are almost identical for LIRGs and 
ULIRGs, explaining the good match between SFRtir 
and SFR R ,o9 at z>2. 

In summary, we conclude that whereas SFRa-ho6 and 
SFRtir are roughly consistent within ^0.3 dex ( mod- 
ulo a constant offset), regardless of the models used to fit 
the IR-SED, the values of SFRbos and SFRr 9 present 
systematic deviations with respect to those that are not 
consistent within the typical rms. Moreover, these dif- 
ferences are not constant, but present a dependence of 
both redshift and SFR. As a result, large systematic off- 
sets (of ±0.5 dex) with respect to SFRtir are expected 
at certain redshifts, e.g., ASFR^+0.5 and -0.5 dex for 
SFRro9 and SFRbos a t z~l-4, respectively. 

5.3.2. Analysis of IR-based SFRs: Best-effort vs. MIPS 

24 firn 

Here we study the impact on the IR-based SFRs of 
modeling the IR-SED with limited photometric data. For 
that matter, we quantify the differences in the SFRs es- 
timated with each of methods compared in the previ- 
ous section using the sample characterized with MIPS 
24 /im data and with 8, 24 and 70 /im data, i.e, the best- 
effort sample (note that <2% of the sample in detected in 



MIPS 70 /*m for ~20% in MIPS 24 /im). In principle, the 
inclusion of additional mid-IR fluxes must improve the 
quality of the estimates given that there is a better sam- 
pling of the IR SED from which better k-corrections to 
the monochromatic luminosities can be obtained. On the 
downside, the spectral range probed by the MIPS bands 
gets narrower with redshift, and the 24 /im channel shifts 
progressively into to PAH region, where models are more 
uncertain and different libraries differ significantly. Also, 
as the observed 70 /im moves further away from the tip 
of the IR-emission (^100/im), the uncertainty in the ex- 
trapolated L(TIR) increases. 

Figure [20] shows the ratio SFR(24) /SFR(8,24,70) as 
a function of redshift for each of the different methods 
to estimate the SFR, except for SFRro9, that only de- 
pends on the observed flux at 24 /im and the redshift) 
. The color code indicates four different bins of infrared 
luminosity. The colored stars with error bars depict the 
median value and 1 er of the redshift and the ratio of 
SFRs(24/8,24,70) for different luminosity bins. 

The upper-left panel of Figure [501 shows that 
SFR T ir(24) and SFR T ir(8,24,70) are ~ler consistent 
within ~0.20dex up to ULIRG luminosities, showing a 
small offset (mostly at z>0.5) in SFRtir (24) towards 
underestimating the SFR by ASFR=-0.05±0.20dex. On 
the other hand, ULIRGs (typically at z>1.5) present val- 
ues of SFRtir (24) larger than SFR T ir(8,24,70) with an 
average difference of ASFR=0.15±0.40dex. This is con- 
sistent with the results found for ULIRGs at this red- 
shift by several authors, who report excesses of a factor 
of 2-10 in the SFR s estimated from MIPS 24 /im only 
(iDaddi et all l2007t iPapovich et~aTI 120071 iRigbv et all 
2008) . Note that estimating the IR- luminosities for these 
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with error bars depict the median value and 1 a of the ratio of SFRs and the redshift per luminosity bin. 



galaxies based on MIPS 24 /jm data alone is intrinsically 
difficult as this band is probing the most variable region 
of the IR-SED, featuring emission from PAHs and Sili- 
cate absorptions. In fact, further motive for these dis- 
crepancies could associated to a change in the relative 
strength of these components in high-z galaxies with re- 
spect to the local templates, particularly for the ULIRG 
templates. 

Figure [21] present further evidence of this issue. The 
left and right panels of the Figure show the rest-frame 
SED normalized to the flux at 24 /im for galaxies at 
z=0.7 and 5=1.5 in different L(TIR) ranges. The red and 
black lines are the dust-emission templates of Rie ke et al.l 
(2009). The black templates are those corresponding 
to the IR-luminosity range shown in the legend. Note 
that we have selected the redshift ranges and template 
normalization with the specific aim of stressing the dif- 
ferences in warm-to-cold dust colors between the R09 



models and the actual observations. The same differ- 
ences apply to other template sets. The vertical line to 
the right of the MIPS 70 /im data depicts the median 
and 1 (7 of the distribution of MIPS 70 fim fluxes (nor- 
malized to 24 /im, i.e., the S70/S24 color). At z^0.7, 
the models for a L(TIR)=10 la5 ~ 110 L nicely predict 
the actual colors observed for galaxies (the templates 
plotted in black match the median and la range of 
observed colors). However, at z^l.5, galaxies present 
smaller colors than what the models for the appropri- 
ate luminosity range predict. This suggest that the ex- 
cess in SFRtir (24) could be related to a difference in 
spectral shapes for ULIRGs at high redshift in compari- 
son with local ULIRGs. Either due to the strength of 
the PAH and the Silicate features (jElbaz et alj l2010t 
iPerez-Gonzalez et al.ll2010T ) or due to additional cont in- 
uum emission by a obscured AGN (jDaddi et al.ll2007t ). 
The lower panel of Figure [2D] shows 
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Fig. 21.— Rest-frame SEDs of galaxies detected in MIPS 24 and 70 /im at 0.6<z<0.8 with 10 10 - 5 <L(TIR)<10 110 (left panel) and 
1.3<z<1.6 with 10 11 - 5 <L(TIR)<10 12 (LIRGS) (right-panel). The fluxes are normalized to the flux in MIPS 24 /im. The solid lines depict 
the templates of Rickc ct al. (2009): black lines are for templates with L(TIR) within the corresponding IR-luminosity range, red lines are 
for the rest of the templates. The vertical line and marks on the right of the MIPS 70 /im fluxes show the median, quartiles and 1 a of the 
distribution of observed fluxes. 

median values of the ratio of SFRs per bin of luminosity 
(and redshift). 

The relative differences with respect to 
SFRtir(8,24,70) remains mostly unchanged with 
respect to what is shown in Figure rTS] for estimates 
based on 24 /mi data. For SFRa-ho6 the most no- 
ticeable differences is that the overall rms increases 
by ^50% for the highest SFRs and the median ra- 
tio for SFR T iR(8,24,7O)>lOOM yr- 1 decreases to 
ASFR=-0.12±0.12dex. For SFR B08 , the median ratio 
SFRbos/SFRtir presents a smaller tilt and a significant 
increment of the rms with respect to the values in the 
MIPS 24 fxm sample. Finally, SFR RO g, which only 
depends on the observed flux in MIPS 24 /mi, presents 
the same trend at low-z as in Figure 1191 However, 
at z^2 it tends to overestimate SFRtir(8,24,70) be- 
cause the former present similar values to SFRtir(24) 
that, as shown in the previous section, overpredicts 
SFR T ir(8, 24,70) at z~2. 

As the relative trends between the estimates of the 
SFR(8,24,70) has barely changed, our conclusions from 
Section l5.3.1l still apply, i.e, the values of SFRa-ho6 are 
roughly consistent with those of SFRtir, but the SFR 
based on the calibrations of B08 and R09 present sys- 
tematic deviations with respect to these that can be sig- 
nificant (up to 0.75 dex) in certain SFR and redshift 
intervals. In addition, we find that the typical rms of 
the comparison of SFR(8,24,70) estimates is 20% larger 
with respect to the previous comparison based on 24 /im 
data. This is not surprising considering that Figure [19] 
shows only functional relation of each method but none 
of the uncertainties attached to the fit of data. 

Note that the comparisons shown in Figure HUl and |2"21 
only illustrates the expected uncertainty budget associ- 
ated to the use of different SFR estimates, but they do 
not demonstrate that any of them provides intrinsically 



SFR(24)/SFR(8,24,70) for the B08 recipe. The ra- 
tio of SFRs presents just the opposite trend of what we 
find for SFRtir, i.e., the offset and rms of the compari- 
son are larger at lower redshifts and almost non-existent 
(< 0.1 dex) at high-z. This is not surprising considering 
that MIPS 24 /im shifts towards 8 /mi with increasing 
redshifts, reducing the impact of the k-corrections. As 
a result, SFRros is nearly insensitive to the inclusion 
of 70 /mi data at z~2. Note however that this does 
not mean that it is a better estimation of the SFR. At 
z<0.5, SFRbo 8 (24) is larger than SFR B08 (8,24,70) by 
ASFR=0.18±0.23dex. 

Finally, the right panel of Figure [20] shows the com- 
parison of SFR(24)/SFR(8,24,70) for the A-H06 recipe. 
The overall trends are analogous to those observed for 
SFRtir but with a larger scatter (~0.30dex), i.e., the 
offset and rms increases with redshift up from ASFR=- 
0.10±0.26dex at z<l to ASFR=0.19±0.47dex at z=2-4. 

5.3.3. Analysis of IR-based SFRs: Best-effort sample 

In this Section we compare again the IR-based SFRs 
obtained with different methods, but this time for the 
best-effort sample, i.e, with estimates based on 8, 24 and 
70 /mi data. We also present a comparison of IR-based 
SFRs to galaxies in common with other authors count- 
ing with better IR SED coverage (e.g., MIPS-160 or IR- 
spectroscopy) and therefore more reliable SFRs. 

Figure [22] shows a comparison of the SFRs(8,24,70) 
obtained with each method with respect to 
SFRtir(8,24,70) (our reference value), as a func- 
tion of SFR (left-panel) and redshift (right-panel). The 
color code is the same as in Figure [12] To simplify 
the comparison to the results of the previous section, 
the Figure also shows the ratio of SFRtir (24) to 
SFRtir(8,24,70) i.e., basically the values shown in the 
upper-left panel of Figure [20] The blue line joins the 
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Fig. 22. — Comparison of different IR-based SFR(8,24,70) indicators with respect to SFR TIR (8,24,70) as a function of SFR TIR (left) and 
redshift (right) for galaxies detected in MIPS 24 /im (f(24 /jm)>60/ijy) and MIPS 70 /an (f(70 /im)>3500/ijy) . The color code is the same 
as in Figure [T9l Here in also show the median values and 1<t of the comparison to SFRbos (g re y stars with error bars). The blue line join s 
the median values of SFRtir(24)/SFRtir(8,24,70) per luminosity bin (similarly to t he colored stars in the to p-left panel of Figure IT9t . 
The black stars de pict the ratio of SFRs for galaxies in common with the sample of Svmconidis ct al. (2008, filled stars at z<1.2) and 
IHuang ct al. (200l|, open stars at z~2). The IR-based SFR for these sources is one of the most accurate available at the moment, and it is 
in overall good agreement with our estimates with SFRtir. 

more accurate results. Nonetheless, having MIPS 70 /mi 
data to constrain the shape of the IR SED, it is reason- 
able to assume that the values of SFRtir(8,24,70) would 
provide more reliable values than the other 3 methods. 
In order to verify this statement and to asses the accu- 
racy of SFRtir (8,24,70), we compare the SFRs to the 
results from other authors based on better photometric 
datasets. In particular, we c ompa re our SF Rs against the 
values of lSymeonidis et all (|2008l S08) and lHuang et al.l 
( 2009). The latter studied the SFRs of a spectroscopic 
sample of high-z (z^l.9) galaxies with strong IR-emission 
(f(24)>0.5 mJy). For these galaxies, the authors provide 
accurate SFRs estimated from a very detailed coverage 
of the IR SED including Spitzer/IRS spectroscopy and 
data at 24, 70 and 160/zm, 1mm and 1.4 GHz. In S08 
the authors describe the IR properties of a 70 /jm se- 
lected sample restricted to galaxies detected at 160 /im 
and having reliable spectroscopic redshifts ranging from 
0.1<z<1.2 (z=0.5). 

The sources in common with S08 and Huan g~et al.l 
(2009) are shown in the right panel of Figure [25] as black 
stars (open and closed, respectively). In addition, we 
show as red stars those sources which were poorly fitted 
to the models of CE01 and DH02 in the work of S08. We 
find that our values of SFRtir (8,24,70) for the majority 
of the z<1.2 galaxies tend to underestimate the SFRs 
of S08 with a median (considering only black stars) dif- 
ference and scatter of ASFR=0.09±q ^ dex. However, 
there is small group of sources for which the SFRs are 
systematically underestimated by ^0.5 dex or more (red 
stars). In S08, the authors showed that, for these galax- 
ies, the IR SED fitting to the models of CE01 and DH02 
severely under-fitted the data at 160/im, whereas the 



models of ISiebenmorgen fc Kriigell (120071 SK07~) allowe d 
a better fit to the data (see also lSvmeonidis et alJ feOlO). 
As a result, the values of L(TIR) obtained from the fit to 
models of CE01 and DH02 would be systematically lower 
than the estimates for SK07. These strong discrepancies 
in the fitting of CE01 and DH 02 models do not seem 
to be the usual scenario (e.g., iKartaltepe et alj [2010). 
although some issues fitting the MIPS 160 /im fluxes of 
local galaxies wit h the models of CE01 has been reported 
dNoll et alJl2009h . 

A possible explanation for this issue could be related 
to fact that S08 makes use of the 4 IRAC bands in the 
SED fitting. As a result, these bands contribute signifi- 
cantly to the x 2 (more than the MIPS bands), whereas 
they only represent a minimum fraction of the total IR 
luminosity. Nonetheless, some intrinsic differences in the 
IR-SED of MIPS 160 /im selected samples are expected, 
given that these are usually biased towards cold galaxies, 
i.e., galaxies with a relatively large (and probably not 
very frequent) cold dust content in comparison with the 
amount and emission of the warm dust featured in the 
models of CE01 and DH02. In summary, the differences 
in the SFR(8,24,70) with respect to the values of S08 for 
these sources are most likely the result of combined SED 
modeling issues and selection effects. 

For the rest of the sources, the ~0.10 dex offset in 
SFR(8,24,70) towards underestimating the values of S08 
is in g ood agreement with the results of lKartaltepe et al] 
(|2010L K10) for a sample of galaxies selected at 70 /mi 
(and counting with 160 /im data for ~20% them). The 
authors indicate that the estimates of L(TIR) for 160 /im 
detected sources computed without fitting that flux can 
be underestimated by up to 0.20 dex at z<l and ~0.3dex 
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at higher redshifts. The authors also point out that these 
effect could be related to a bias in 160/im selected sam- 
ples towards selecting intrinsically cooler objects (as op- 
posed to 70 /im selections). A similar statement is made 
in S08 based on the f 60/70 flux rat ios of their sample. 

The comparison to the SFRs of iHuang etafl (|2009| ) 
for galaxies at z— 2 present too few sources to provide 
a significant result. However, the overall results are in 
relatively good agreement within 0.3 dex. We find again 
that SFR T ir(8,24,70) slightly under-predicts the SFRs 
of some galaxies, consistently with the conclusions of 
K10 for 160/im detected galaxies. Moreover , we find that 
the va lues of SFRtir(24) for the galaxies of lHuang et al.l 
( 2009) overestimates SFR(8,24,70) by a factor of ~4, fol- 
lowing the trend shown by the blue line. 

Finally, we further check the accuracy of our estimates 
at high-z, comparing them to the values of SFRtir(24) 
correc ted with the empirical relation of Papovic h et al.l 
(2006). This correction was conceived to mitigate the ex- 
cess in the IR-SFRs of high-redshift galaxies estimated 
from 24 /im data. The correction was computed by 
matching the SFR(24) to the SFRs estimated from the 
average stacked fluxes in MIPS 24, 70 and 160 fan of a 
sample of z— 2 galaxies. The green line in Figure [2"2l joins 
the median ratios of SFR T ir,corr(24)/SFR T ir(8,24,70) 
as a function of redshift. The overall results are 
that the values of SFRtir,corr(24) are a factor of 
—3-4 lower than SFRtir (24) at z— 2. As a result 
these estimates are also slightly lower than our pre- 
dictions for SFR T ir(8,24,70). Nonetheless, the val- 
ues of SFRtir,corr(24) also agree within 0.3 dex with 
the SF Rs of the galaxies in common with IHuang et al.l 
([2009T) . 

5.4. Summary of the SFRs 

The accuracy of the SFRs estimated from IR trac- 
ers up to intermediate redshifts has been demonstrated 
by the good agreement with the estimates based on 
oth er tracers such us dust corrected UV/optical indica- 
tor (jlglesias-Paramo et alj2007t ISalim et al.l2007l[2009t ). 
On the other hand, the systematics effects in the IR- 
based SFRs of the most luminous galaxies (ULIRGS) at 
high redshift are quite significant. Some of these issues 
arise from the assumptions made in the estimation of IR- 
based SFRs, such as the validity of the local templates 
at high redshift or the contribution of obscured AGNs 
to the IR luminosity. However, the most relevant issues 
arise from the lack of enough data to constrain the full 
IR SED. Particularly for studies based on 24 /im data 
alone. Nonetheless, the breadth and quality of the MIPS 
24 /im data ensures that it will continue leading multiple 
studies of IR-based SFR for the foreseeable future. Thus, 
quantifying the systematic effects between the SFRs(24) 
computed with different methods, and the differences in 
the SFRs(24) with respect to the SFRs computed from 
more IR data, provides a useful information. 

Our analysis shows that, although the values of 
SFR T ir(24) are consistent with those of SFR A -ho6(24) 
within 0.3 dex (the usual uncertainty quoted for IR-based 
SFRs) the values of SFR B08 (24) and SFR R09 (24) can be 
significantly deviated (up to ±0.5 dex) with respect to 
SFRtir(24) for certain redshift and luminosity ranges. 
The differences in the SFRs obtained with these methods 
remain mostly unchanged for SFR(8,24,70), and we find 



that the discrepancies in SFRbos or SFRro9 with respect 
to SFRtir(24) do not provide a better agreement to the 
SFRs of other authors computed from very detailed IR 
photometric data. Therefore, out of the four methods 
to estimate the IR-SFR discussed here, SFRtir present 
(after accounting for intrinsic systematics) the more ac- 
curate results. 

From the analysis of sample of MIPS 70 /im detected 
galaxies, we find that SFRtir(24) is reasonably consis- 
tent with the values of SFR T ir(8,24,70) up to ULIRG lu- 
minosities (typically at z<1.4) showing only a small devi- 
ation towards underestimating SFR(8,24,70) by 0.05 dex 
with an rms of 0.2 dex. However, at z>1.5 the agree- 
ment is significantly worse. The values of SFRtir (24) 
tend to overestimate SFR(8,24,70) by a median value of 
0.15±0.40 dex. As already pointed out by other authors, 
the best approach to solve this issue is to apply a correc- 
tion factor that redu ces the estimated valu es at high-z 
(|Papovich et al.ll2006l or lSantini et al.ll2009t ). 

The comparison of SFRtir (8,24,70) (our best-effort 
SFRs) to the SFRs computed by other authors based 
on a better IR photometric coverage (including MIPS 
160 /xm) also shows an excellent agreement, proving that 
these estimates are robust. The overall results are consis- 
tent within 0.3 dex presenting only a small systematic de- 
viation in SFRtir(8,24,70) towards underestimating the 
values including MIPS 160 /mi data by -0.09 dex (mostly 
z<1.2 galaxies). Note that since this comparison is re- 
stricted to MIPS 160 /im detected sources the could be 
some selection effects, and thus this offset might not ap- 
ply for all galaxies (see e.g., the results of K10 based 
on stacked fluxes in MIPS 160 /im for a 70 /im selected 
sample). 

6. DATA ACCESS 

All the data products for the 76,936 IRAC 3.6+4.5 fim- 
selected ([3.6]<23.75) sources in the EGS are presented 
here. These include: (1) the photometric redshift cata- 
log containing the estimates with Rainbow, EAZY and 
from 106, when available (Table [7]); (2) the stellar mass 
catalog containing the values estimated with each of the 
different modeling configurations described in § U (Ta- 
ble [SJ); (3) the SFR catalog containing the UV- and IR- 
based SFRs obtained with the different methods and cal- 
ibrations discussed in §[5] (Table |9]). A table containing 
the UV-to-FIR SEDs for all these sources is presented in 
Paper I. The number of objects and unique identifier of 
this table and the tables presented in the following is the 
same. 

A larger version of these catalogs containing all galax- 
ies down [3.6]<24.75 (3ct limiting magnitude) are avail- 
able through the web utility Rainbow Navigator 2 (see 
paper I for a more detailed description), that provides 
a query interface to the database containing all the data 
products of the multiple Rainbow tasks that we have used 
in the papers. Rainbow Navigator has been conceived to 
serve as a permanent repository for future versions of the 
data products in EGS, and also to similar results in other 
cosmological fields (such as GOODS-N and GOODS-S, 
presented in PG08). 

6.1. Ta6/e[?|' Photometric redshift catalog 
2 http://rainbowx.fis.ucm.es 
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These are the fields included in Table [7] 

• Object: Unique object identifier starting with 
iracOOOOOl. Objects labeled with an underscore 
plus a number (e.g, irac000356_l) are those iden- 
tified as a single source in the IRAC catalog built 
with SExtractor, but deblended during the photo- 
metric measurement carried out with the Rainbow 
software (see § 12. ip . Note that, although the cata- 
log contains 76,936 elements, the identifiers do not 
follow the sequence iracOOOOOl to irac076185. This 
is because the catalog is extracted from a larger ref- 
erence set by imposing coordinate and magnitude 
constraints. The table is sorted according to this 
unique identifier. 

• a, 5: J2000.0 right ascension and declination in de- 
grees. 

• zphot-peak: Maximum likelihood photometric red- 
shift. 

• zphot-best: Probability weighted mean photometric 
redshift. This is the value of z p hot used along the 
paper. 

• zphot-err. ler uncertainty in the photometric red- 
shift as estimated from the zPDF. 

• zphot-EAZY: Photo metric redshift estima ted us- 
ing the EAZY code (|Brammer et al.ll2008f) on our 
SEDs with the default templates and including the 
JT-band luminosity prior. 

• Qz: Estimate of the quality of the photometric red- 
shifts computed with EAZY. Relia ble photometric 
redsh ifts present values of Qz<l (|Brammer et al.l 
[2001 . 

• zphot-1 06: Photometric redshift from lllbert et al.l 
(|2006a[ ). These are only available for galaxies in 
the main region. 

• zspec: Spectroscopic redshift (set to -1 if not avail- 
able). 

• qflag: Spectroscopic redshift quality flag from 1 to 
4. Sources with qflag>3 have a redshift reliability 
larger than 80%. 

• N (bands): Number of photometric bands used to 
derive the photometric redshift. 

• Stellarity: Total number of stellarity criteria satis- 
fied. A source is classified as a star if it satisfies 3 
or more criteria. A description of all the stellarity 
criteria and the accuracy of the method is given in 
§ 5.4 of Paper I. 

6.2. Tabled stellar mass catalog 

The stellar masses are estimated from the same tem- 
plates used to compute the photometric redshifts. These 
templates were computed using several combination of 
stellar population synthesis library, IMFs and extinc- 
tion laws. Our reference stellar masses are those ob- 
tained with [P01,SALP,CAL01] (see jj3J}. We provide 



2 different stellar mass estimates based on these tem- 
plates depending on the redshift used during the fit- 
ting procedure, namely: zphot-best and zspec. In ad- 
dition, we additional estima tes obtained with: 1) the 
stellar population mode ls of iBruzual fc Charlotl (2003, 
BC03). PMarastonl (|200H M05) and Cha r iot fc Bruzual 
(200 9; CB09) 2) Th e IMFs of iKroupal (|200ll KROU) 
and IChabrierl (120031. CHA B), and 3) the dust extinc- 
tion law~o?lC^rlQr&M3 (l2000l CF00). In Tabled 
we give 6 additional stellar mass estimates obtained un- 
der different stellar population modeling assumptions, 
namely: [P01,KROU], [BC03,CHAB], [M05,KROU], 
[CB09,CHAB], [CB09,SALP] and [P01,SALP,CF00]. 
The extinction law in all cases except for the last is 
CAL01. For these stellar mass estimates we use z-fit, 
which is equal to zphot-best unless zspec is available. 
These are the fields included in Table [5] 

• Object: Unique object identifier as in the photo- 
metric catalog. 

• a, 8: J2000.0 right ascension and declination in de- 
grees. 

• Mass(best): Stellar mass [log M©] with the 
associated uncertainty estimated with zphot- 
best using our default modeling parameters 
[P01,SALP,CAL01]. 

• Mass(zspec): Stellar mass [log M Q ] with the asso- 
ciated uncertainty, estimated with zspec using our 
default modeling parameters [P01,SALP,CAL01]. 

• z-fit: Value of the photometric redshift used 
during the SED fitting with the [P01,KROU], 
[BC03,CHAB], [M05,KROU], [CB09,CHAB], 
[CB09,SALP], [P01,SALP,CF00] models. It is 
equal to zphot-best unless zspec is available. 

• Mass(P01,KROU): Stellar mass [log M Q ] with the 
associated uncertainty, estimated with the model- 
ing parameters [P01,KROU,CAL01] and zphot-fit. 

• Mass(BC03,CHAB): Stellar mass [log M Q ] with the 
associated uncertainty, estimated with the model- 
ing parameters [BC03,CHAB,CAL01] and zphot- 
fit. 

• Mass(M05,KROU): Stellar mass [log M Q ] with the 
associated uncertainty, estimated with the model- 
ing parameters [M05,KROU,CAL01] and zphot-fit. 

• Mass(CB09,CHAB): Stellar mass [log M ] with the 
associated uncertainty, estimated with the model- 
ing parameters [CB09,CHAB,CAL01] and zphot- 
fit. 

• Mass(CB09,SALP): Stellar mass [log M Q ] with the 
associated uncertainty, estimated with the model- 
ing parameters [CB09,SALP,CAL01] and zphot-fit. 

• Mass(P01,CF00): Stellar mass [log M ] with the 
associated uncertainty, estimated with the model- 
ing parameters [P01,SALP,CF00] and zphot-fit. 
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6.3. Tabled SFR catalog 

The unobscured UV-SFR is obtained from 
the best-fitting optical template modeled with 
[P01,SALP,CAL01]. The rest-frame IR luminosi- 
ties and IR-based SFRs are computed either from the 
average value of the best fitting temp l ates f rom the dust 
emission mod e ls of [Ch ary fe Elbazf (|200l CE01) and 
iDale fe Heloul ([20021 DH02) to MIPS 24 /im data only, 
or the average value of the best fitting dust emission 
models of CE01, DH02 and iRieke efalT (|2009l R09) to 
IRAC-8.0, MIPS 24 and 70 /im data. Only fluxes at 
rest-frame A >5 /im are considered in this method. In 
both cases we use zphot-fit, which is equal to zphot-best 
unless zspec is available. For sources un-detected in 
MIPS 70 (im at z>0.6 both methods provide similar 
results modulo the effect of the R09 templates. 

Note that the SFRs has been computed for all the 
MIPS 24 /jm and 70 /jm detections, but only sources with 
f(24)>60 /iJy and f(70)>3.5mJy (the 5a detection limit) 
are discussed in § [Sj In addition, sources un-detected 
in MIPS 24 /im are fitted using an upper limit value of 
f(24)=60 /iJy. In this cases the quoted L(TIR) and SFRs 
are negative values. 

• Object: Unique object identifier as in the photo- 
metric catalog. 

• a, 5: J2000.0 right ascension and declination in de- 
grees. 

• f(24),f(70): Observed flux [/iJy] and uncertainties 
in MIPS 24 and 70 /im. 

• z-fit Value of the redshift used during the IR SED 
fitting. It is equal to zphot-best unless zspec is avail- 
able. 

• L(TIR,24): Total IR luminosity [log Moyr^ 1 ], cal- 
culated by integrating the (average) dust emission 
model from 8 /im to 1000 /im . This value is com- 
puted by fitting the observed flux in MIPS 24 /im 
to the models of CE01&DH02. 

• SFR . 28 : Unobscured UV-based SFR [Moyr" 1 ] es- 
timated from the rest- frame luminosity at 0.28 /im 
interpolated in the best-fit optical template, 
i/Ly (0.28), using the Kennicutt (1998) calibration. 

• SFR T ir(24): IR-based SFR [Moyr^ 1 ], estimated 
from L(TIR) using the calibration of Kennicutt 
(1998). This value is computed by fitting the 
observed flux in MIPS 24 /im to the models of 
CE01&DH02. 

• SFR CE0 i(24): Same as SFR T ir(24) but fitting the 
MIPS 24 /on data to the models of CE01 only. 

• SFR B os(24): IR SFR [Mayr" 1 ] estimated from the 
rest-frame monochro matic luminosity at 8a m us- 
ing the calibration of iBavouzet et ail (|2008! ). This 
value is computed by fitting the flux in MIPS 24 /im 
to the models of CE01&DH02. 

• SFR A -ho6(24): Total SFR [Moyr" 1 ] esti- 
mated from the rest-frame monochromatic 
luminosity at 24 fim using the calibration of 



lAlonso-Herrero et al.l (|2006l ). This value is com- 
puted by fitting the observed flux in MIPS 24 /im 
to the models of CE01&DH02. Note that to 
obtain the IR-SFR part of this value, the unob- 
scured UV-SFR must be subtracted according to 
Equation 01 

• SFR R09 (24): IR SFR [Moyr^ 1 ] estimated from the 
observed flux i n MIPS 24 /im and the redshift using 
the formula of IRieke et all (|2009l eq 14). 

• L(TIR) ,SFRtir,SFR C eoi ,SFR B o8,SFR A -ho6 (8,24,70) . 
Same as the previous values but fitting the IR 
SED with IRAC-8.0, MIPS 24 and 70 /im data 

to the models of CE01, DH02 and R09. Note 
that SFRro9 has been omitted because its value 
is independent of the flux in IRAC-8.0 nor MIPS 
70 /im. 

7. SUMMARY 

In this paper, and the companion (Barro et al. 2010a; 
Paper I), we have presented an IRAC-3.6+4.5/im selected 
sample in the Extended Groth Strip characterized with 
UV-to-FIR SEDs. The photometric catalog includes the 
following bands: far-UV and near-UV from GALEX, 
u*g'r'i'z' from the CFHTLS, u'gRiz from MMT and 
Subaru observations, BRI from CFHT12k, V606, *8i4, 
Jno and -Hi6o from HST , JK from Palomar Observa- 
tory, CAHA and Subaru data, and [3.6]-to-[8.0], 24 /im, 
and 70 /im data from Spitzer IRAC and MIPS surveys. 
Our catalog contains 76,936 sources down to a 85% com- 
pleteness level ([3.6]<23.75) over an area of 0.48deg 2 . In 
addition, we have cross-correlated our sample with the 
redshift catalog from DEEP2, and with X-ray and VLA- 
20cm radio data. 

Paper I presented the data, the procedure to measure 
consistent UV-to-FIR photometry using our own dedi- 
cated software (Rainbow), and the analysis of the multi- 
band properties of the sample. We showed that the SEDs 
present the level of consistency required to characterize 
the intrinsic stellar populations of the galaxy. In this pa- 
per, we have presented a galaxy-by-galaxy fitting of the 
UV-to-FIR SEDs to stellar population and dust emission 
models. From the best fitting optical and IR templates, 
we have estimated: (1) photometric redshifts, (2) stellar 
masses, and (3) SFRs. Then, we have analyzed in detail 
their accuracy and reliability with respect to different 
parameters. In the following we present the summary of 
the most important results of this analysis, organized by 
parameter. 

Photometric redshifts (z p hot) were estimated from the 
comparison of the UV-to-NIR SEDs to stellar popula- 
tion and AGN templates. This comparison was car- 
ried out with our own dedicated software (within the 
Rainbow package) using y 2 minimization algorithm (see 
iPerez-Gonzalez eta l. 2008 for more d etails), and with 
the EAZY code (jBrammer et al.l 120081 ). These are our 
main results about photometric redshifts: 

• Two new features have been included in the Rain- 
bow photometri c redshift code over the previ ous im- 
plementation in lPerez-Gonzalez et al.l (|2008j ) to im- 
prove the quality of the estimates: (1) a zero-point 
re-calibration of the observed photometry, and (2) 
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the use of template error function as a weight term 
in the SED fitting procedure. Both features are 
computed simultaneously and iteratively based on 
the comparison of observed and synthetic photom- 
etry in a spectroscopic control sample. The results 
show an overall good agreement between observa- 
tions and templates. The zero-point corrections are 
typically <0.1 mag and converge after a few itera- 
tions. The overall rms in the residual is a factor of 
^2 the median photometric uncertainty. The most 
noticeable discrepancies present at A>3 /im and (to 
a lesser extent) around 250nm. These are the result 
of limitations in the stellar templates in the NIR 
range, and a possible excess in the s trength of the 
dust a ttenuation with respect to a Ca lzetti et al.l 
(2000) extinction law, respectively. 

The comparison of our photometric redshifts 

to 7636 secu re spectrosco pic redshifts from 

DEEP2 and iSteidel et al.l (l200l LBGs at 
z>3) shows an overall accuracy of <7ncmad= 
Az-median(A2 



1.48 x median i 



=0.034 (where 



Az=z pho t - z spcc ) and cr NC MAD=0.046, with rj=2% 
and 3% catastrophic outliers {q defined as the 
fraction of galaxies presenting oncmad>0.2) in 
the EGS main region (covered by the CFHTLS) 
and flanking regions (covered with fewer and 
shallower bands), respectively. The overall scat- 
ter in Zphot does not depend strongly on the 
redshift, presenting a minimum value around 
z=0.5-l ((Tncmad=0.028 and 0.040 in the main 
and flanking regions, respectively), and increasing 
by a factor of ~1.3 at lower and higher redshifts 
(up to z<1.5). At z^3, the z p hot accuracy for 
91 LBGs with secure spectroscopy is reduced to 
cncmad=0.063 and 77=10%. 

• The accuracy of the z p hot is mostly independent 
of the [3.6] magnitude. However, it decreases 
with the optical magnitude from <tncmad=0.030 at 
R=22 to ctncmad=0.060 at R=25. Approximately 
50% of the catastrophic outliers have R>23.5 and 
log(% 2 )>0.6. Approximately 60% of the sources 
with significantly different values of Zb es t and z pca k 
(|^best-^pGak|/(l+z)>0.2) are catastrophic outliers. 

• The Zphot statistics for the 1995 and 262 spec- 
troscopic galaxies detected in MIPS 24 /im 
(f(24)>60 /jjy) and 70 /im (f(70)>3.5 mJy) in the 
main region are similar to the rest of the sample 
with ctncmad =0.033, r)=3%, and er NC MAD=0.045, 
77=1%, respectively. The accuracy for the 142 X- 
ray sources is similar (oncmad=0.038) although 
with larger fraction of outliers (77 =10%), proba- 
bly as a result of some degree of contamination 
by the AGN, for which reliable z p hot are diffi- 
cult to estimate based on stellar templates. The 
worst results are found for a very few (12) power- 
law galaxies (PLGs, identified as obscured AGNs): 
cncmad=0-052 (77 =17%). We also note that 
sources with increasing fluxes in the IRAC bands 
(f[3.6]<f[4.5]<f[5.8]<f[8.o]) makes up for up to 15% of 
the total number of outliers. 



• The z p hot Rainbow are in good agreement 
with those from t he i'-band selected catalog of 
(jllbert et al.ll2006al 106), which overlaps with our 
sample in the main region. For the 5454 galaxies in 
common between the two catalogs with [3.6]<23.75 
and i' <24.5, the accuracy of the z p hot at z<l is 
roughly the same, <tncmad=0.035. At higher red- 
shifts, our larger band coverage (mostly in the NIR) 
provides more accurate results and less severe sys- 
tematic errors and uncertainties. In particular, for 
galaxies at z^3 (the LBG sub-sample), the outlier 
fraction in 106 is 46% for only 9% in Rainbow. Our 
Zphot catalog and the one presented in 106 are com- 
plementary: whereas the NIR-sclccted sample de- 
tected more galaxies at high-z, which are too faint 
in the optical to be included in the 106 catalog, 
the IRAC catalog misses a population of low-mass 
galaxies at z<l which are recovered by the i'-band 
selection in 106. 

• We showed that the photometric catalog provides 
robust SEDs by obtaining a different realization of 
the z P ho t catalog with similar quality using the code 
EAZY (|Brammer et al.ll2008T) . In particular, these 
alternative photometric redshifts are slightly more 
accurate for the sources in the main region, par- 
ticularly at z<0.5, whereas they present a larger 
scatter in the flanking regions. Moreover, these 
Zphot exhibit a slightly larger systematic deviation 
(Az/(l+z)=0.019 and 0.027 in the main and flank- 
ing regions, respectively) than the z p hot computed 
with Rainbow. 

• We further tested the accuracy of our z p hot by 
checking the number densities and z p hot distri- 
butions of a sub-sample of (NIR selected) s-BzK 
(p =5.0 arcmin 2 ; z =1.89), p-BzK (p =0.5 arcmin 2 ; 
I =1.85) and DRG (p =1.4 arcmin 2 ; z =2.47). 
These are in relatively good agreement with the re- 
sults from the literature down to ifyEGA<21. The 
most significant difference is an excess of ~1.5 in 
the density s-BzKs, which could be caused by an 
overdensity of galaxies at z^l.5. 

• The median redshift of the ([3.6]<23.75) sample, 
2=1.2, is consi stent with that of the flux lim- 
ited sample s of fPerez-Gonzalez et al.l (|2008l) and 
lllbert et al.l ((20091) in different fields. 

Stellar masses for the whole sample were obtained in 
a galaxy-by-galaxy basis by fitting the optical-to-NIR 
SEDs to stellar population synthesis models. In addition, 
we analyzed the effects of the choice of different stellar 
population synthesis (SPS) libraries, IMFs and dust ex- 
tinction laws on our estimations. For that, we considered 
a reference set of assumptions to which several combina- 
tions of input parameters were compared. This refer- 
ence stellar masses were obtained with the PEGASE 2.0 
(jFioc fc Rocca-Volmera nsc 1997) stellar population syn- 
the sis models (P01 ) a Sa lpeter (1955) IMF (SALP) and 
the ICalzetti et ail (|2000l CALZ01) extinction law. We 
compared these estimations with t hose obtained with: 
1) th e stellar population model s of iBruzual fc Charlotl 
2001 BC031. PMalastonl (120051 M05) and Cha r iot fc 
Bruzual (2009; CB09); (2) the IMFs of lKroupal (l200l 
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KROU) and Chabricr (20031 CHAB); (3) the dust ex- 
tinction law of ICharlot fc Falll (|2000L CF00). These are 
our main results about stellar masses: 

• From the comparison of the stellar masses esti- 
mated with photometric and spectroscopic red- 
shifts we find a ler uncertainty of —0.2 dex. The 
distribution of stellar masses as a function of red- 
shift for our default modeling assumptions shows 
that 90% of the galaxies present log(M)>10M Q 
at z>2, at the limiting magnitude of our sample 
Q3.6]<23.75). 

• We quantified the impact of the choice of differ- 
ent IMFs in the estimated stellar masses. For the 
CB09 models, we found that the use of a SALP, 
KROU, or CHAB IMF introduces constant off- 
sets (with a very small scatter) in the estimated 
stellar masses: A log(M)[SALP— KROU]=0.19dex 
and Alog(M)[CHAB-KROU]=-0.04dex. For 
the models of P01, the difference for a SALP 
and KROU IMFs depends on the mass, rang- 
ing from Alog(M)=0.03dex for masses lower than 
log(M)=10A4© to 0.13dex above that threshold. 

• We quantified the impact of using different 
SPS codes in the estimated stellar masses. 
We found that the new CB09 models pre- 
dict slightly lower masses than the older ver- 
sion, BC03, by Alog(M)=0.04 ±g;f|dex. Our 
stellar masses estimated with the P01 mod- 
els are on average larger than those obtained 
with the CB09 models (for a KROU IMF) by 
Alog(M)=0.15±0.26dex. The estimates with the 
P01 library are also larger than those with the 
M05 SPS by Alog(M)=0.39±0.34dex. We found 
slightly lower values of this offset for galaxies with 
log(M)>10X Q (-0.30 dex). Our default modeling 
assumptions, [P01,SALP,CAL01], predicts com- 
paratively the largest stellar masses. Accounting 
for all systematic offsets, all models are roughly 
consistent within a factor of 2-3. 

• We quantified the effect of using different treat- 
ments of the dust extinction by comparing the stel- 
lar masses estimated with a CAL01 and CF00 ex- 
tinction laws. The median result is a small system- 
atic deviation of 0.03 dex towards smaller values 
when using CAL01, and a rms of —0.20 dex. This 
suggest that the different treatments of the dust at- 
tenuation do not play a major role in the estimate 
of the stellar masses. 

• The comparison of our results with several stellar 
mass catalogs already published in EGS revealed a 
good agreement despite the differences in the mod- 
eling technique and in the photometric dataset. We 
found a median offset and scatter of Alog(M)=- 
0.07±0.21dex and A log(M)=0.10±0.25dex with 
res pect to the catalogs of st ellar masses publishe d 
bv iBundv et~ai1 (2006) and iTruiilloeFahl (20071 ), 
respectively. 

SFRs were estimated for all galaxies in our sample fol- 
lowing a variety of procedures. First, we calculated the 



unobscured SFR (the star formation which is directly 
observable in the UV/optical) from the observed lumi- 
nosity at 280 nm (SFRyy^s). To get the total SFR of a 
galaxy, the former value must be added to the SFR which 
is not directly measurable in the UV/optical because of 
the extinction by dust. We calculated this SFR from 
the IR data taken by Spitzer/IRAC and MIPS at 24 /xm 
and 70 /mi(if available). The general procedure consist 
on fitting the IR photometry at rest-frame wavelengths 
A >5/xm (usually involving 8, 24 and 70 iim data ) to the 
dust emission templates of iCharv fc Elbazl (f200lL CE01) 
and lDale fc Heloul (20021 DH02) and lRieke et all (20091) , 
but we also performed some test by fitting only MIPS 
24 iim data to the models of CE01&DH02. 

From the best fit to the models, the IR-based SFR for 
each object was estimated with 4 different methods: (1) 
the total infrared luminosity, L(TIR), integrated from 8 
to 1000 [im transformed to a SFR with the factor pub- 
lished bv lKennicuttl (1998T) . (2) the rest-frame monochro- 
matic luminosity at 8 /im (SFRbos) transformed to 
L(TIR) and SFR using the empirical relation described in 
iBavouzet et alj (20081 ) and the Kennicutt factor ( 3) the 
empirical relation given in lAlonso-Herrero et al.1 (2006ft 
between the rest-frame monochromatic luminosity at 
24 ^m and the SFR ( SFR A -ho6)- (4) using equation (14) 
of iRieke et all (2009 ), that relates the SFR (SFR R09 ) to 
the observed flux in the MIPS 24 [im band and the red- 
shift. The monochromatic and integrated luminosities 
were computed from the average value of the best fit 
templates. These are our main results about SFRs: 

• We quantified the differences in the IR-based SFRs 
obtained with the four methods based on the fit 
to MIPS 24 ^m data only. The SFRs estimates 
with the models of CE01 and DH02 are com- 
patible within a factor of 2, presenting a max- 
imum difference around z— 1.5. The estimates 
of SFRtir and SFRa_ho6 ar e roughly consistent 

(ASFR 0.18±0.05dex) when the contribution of 

SFRyv^o&s is small. SFRbos gives systematically 
lower values than SFRtir for SFR>20 M yr _1 
and z>l, and higher values for lower redshifts and 
SFRs. The difference exceeds a factor of 5 for 
SFR>1000M o yr _1 . The overall agreement be- 
tween SFRtir and SFRro9 is rather poor, ex- 
cept at z>1.8 where the differences are lower than 
0.05 dex. The reasons for these discrepancies can 
be found in the differences in the relative emission 
of the cold and warm dust, and in the strength of 
the PAH and silicate absorption. These character- 
istics can vary by up to a factor of —5 from one set 
of templates to the other. 

• For each of the methods to estimate the SFR, 
we studied the effect having a better constrained 
IR SED comparing the SFRs computed from 
IRAC+MIPS, SFR(8,24,70), and just MIPS 24 fim, 
SFR(24). At low-z, the median values of 
SFRtir (24) and SFR A -Hoe(24) tend to underes- 
timate SFR(8,24,70) by 0.05 and 0.10 dex, respec- 
tively, with an rms of —0.2-0.3 dex. At z— 2, the 
estimates from these two methods based on 24 /im 
data only are on average —0.20 dex larger than 
those obtained with SFR(8,24,70). SFR B08 (24) 
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presents the opposite trend, giving larger values 
than SFRbos (8,24,70) at z<l (up to 0.18 dex), but 
remaining mostly unchanged at higher redshifts. 

• The relative differences between each of the meth- 
ods to estimate the IR-based SFRs with respect 
SFRtir (described in the first item) remains 
mostly unchanged when using SFRs(8,24,70). 
The values of SFR T ir(8,24,70) (best-effort esti- 
mate) for a sam p le of galax ies in common wit h 
iSvmeonidis et all (|2008f ) and iHuang et al.l (|2009|) . 
who counted with more photometric fluxes in the 
FIR and (sub-)mm range, presented a relatively 
good agreement within ~0.3dex. At z<1.2 we find 
a small deviation of 0.09 dex in SFRtir (8,24,70) 
towards underestimating the SFR of MIPS-160 de- 
tected galaxies. Larger discrepancies, up to 0.5 dex, 
might arise for individual galaxies due to the use 
of different template sets. 

In the context of studies of galaxy evolution, our cata- 
log provides a self consistent sample with a very detailed 
characterization of the systematic uncertainties suitable 
for multiple scientific purposes. It is also an alternative 
to other catalogs providing only photometry, redshifts or 
stellar parameters alone. Furthermore, our photometric 
catalog itself provides a reference point for independent 
analysis of the stellar populations. 

The multi-band photometric catalog presented in Pa- 
per I, jointly with the photometric redshifts and es- 
timated stellar parameters presented here are publicly 
available. We have developed a web-interface, named 
Rainbow Navigator 3 , that provides full access to the 
imaging data and estimated parameters and allows sev- 

3 http://rainbowx.fis.ucm.es 



eral other data handling functionalities. 
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TABLE 7 

The IRAC-3.6+4.5 /jm sample: Photometric redshifts 



Object 


a 


<5 


zphot-best 


zphot-EAZY 


zphot-106 


zspec 


qflag 


zphot-err 


Qz 


N(band) 


Stellarity 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


irac003270.1 


215.43910540 


53.08468920 


1.09 


1.06 


1.04 


0.00000 


2 


0.07 


0.31 


16 





irac003278 


215.42614011 


53.09447161 


0.49 


0.39 


0.58 


0.00000 





0.20 


2.58 


16 


7 


irac003291_l 


215.44058360 


53.08123980 


0.83 


0.80 


0.84 


0.85700 


4 


0.02 


0.11 


16 





irac003310 


215.42129738 


53.09430607 


0.08 


0.14 


0.15 


0.00000 





0.07 


0.13 


18 





irac003313 


215.43553774 


53.08200958 


1.06 


1.01 


1.09 


0.00000 





0.08 


0.37 


16 


1 



Note. — (1) Object unique identifier in the catalog. 
(2,3) Right Ascension and Declination (J2000) in degrees. 

(4) Probability weighted photometric rcdshift. This is our default value of photo metric redshift for SED fitting based estimates. 

(5) Photometric rcdshift estimated with the code EAZY ([Brammcr ct al. 200$) using the default template configuration and the A"-band luminosity prior applied to the 
[3.6] band. The input photometric catalog is the same as for t he other redshifts. 

(6) Photometric redshifts as estimated in jllbcrt et al.j < j2006a|) from the (5-band) % selected catalog of the CFHTLS. This catalog overlaps with the IRAC sample in the 
central portion of the mosaic (52.16° < S <53.20° & 214.04° < a <215.74°). 

(7) Spectroscopic rcdshift determination drawn from DEEP2 dDavis et al.|[2007l ; ~8,000 galaxies) and dSteidel et al .112003 ; LBGs at z>3). 

(8) Quality flag of the spectroscopic rcdshift (4— >99.5%, 3— >90%, 2— uncertain, 1— bad quality). Only redshifts with qflag>2 have been used in the analysis. 

(9) Uncertainty in zphot-bcst(4) estimated from the 1 a width of the probability distribution function . 

(10) Reliability parameter of the photometric redshift estimated with EAZY (sec Brammcr ct al. 2008 for more details); Good quality redshifts are in general Q z <1. 

(11) Number of different photometric bands used in to estimate the photometric rcdshift with Rainbow, column (4). 

(12) Sum of all the stellarity criteria satisfied (sec §5.4 of Paper I). A source is classified as star for Stcllarity>2. 

(This table is available in its entirety in a machine-readable in the online version. A portion is shown here for guidance.) 



TABLE 8 

The IRAC-3.6+4.5 /um sample: Stellar Mass estimates 



Object 


a 


5 


M(bcst) 


M(zspec) 


z-fit 


M(P01,KROU) 


M(BC03,CHAB) 


M(M05,KROU) 


M(CB09,CHAB) 


M(CB09,SALP) 


M(P01,CF00) 








M-err 


M-err 




M-crr 


M-crr 


M-crr 


M-err 


M-crr 


M-err 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 








(13) 


(14) 




(15) 


(16) 


(17) 


(18) 


(19) 


(20) 


irac003270.1 


215.43892696 


53.08455063 


9.79 


9.79 


1.09 


9.67 


9.64 


9.61 


9.61 


9.84 


9.98 








0.05 


0.05 




0.04 


0.08 


0.07 


0.08 


0.08 


0.07 


11 acUUOZ i o 


91 ^ <19fi1 401 1 


ccq 0Q/M71 fi1 


1 1 fiQ 


1 1 RQ 




11^ 


1 l.UU 


1 1 ^ 


11 48 


11 74 


1 l 7H 








0.15 


0.15 




0.14 


0.17 


0.16 


0.18 


0.15 


0.17 


irac003291_l 


215.44043562 


53.08128671 


10.84 


10.84 


0.86 


10.97 


10.64 


10.30 


10.57 


10.63 


10.48 








0.03 


0.04 




0.04 


0.04 


0.03 


0.03 


0.04 


0.04 


irac003310 


215.42129738 


53.09430607 


8.80 


8.80 


0.08 


8.61 


8.92 


8.64 


8.79 


9.04 


8.75 








0.06 


0.06 




0.06 


0.08 


0.08 


0.05 


0.05 


0.06 


irac003313 


215.43553774 


53.08200958 


9.94 


9.94 


1.06 


9.83 


9.64 


9.19 


9.40 


9.63 


9.59 








0.04 


0.04 




0.05 


0.04 


0.05 


0.07 


0.07 


0.04 



Note. — (1) Object unique identifier in the catalog. 
(2,3) Right Ascension and Declination (J2000) in degrees. 

(4-13) Stellar mass [log(Mo)] with the associated uncertainty, estimated with our default modeling parameters, [P01,SALP,CAL01] , and zphot-best. 

(5-14) Stellar mass [log(MQ)] with the associated uncertainty, estimated with our default modeling parameters, [P01,SALP,CAL01] , forcing the photometric rcdshift to the spectroscopic value, 
when available. 

(6) Rcdshift used in the fitting procedure of (8,9,10,11 and 12). This redshift refer to zphot-best unless the spectroscopic rcdshift is available; in that case the rcdshift is forced to the spectroscopic 
value. 

(7-15) Stellar mass [log(Mo)] with the associated uncertainty, estimated with the modeling parameters, [P01,KROU,CAL01], and z-fit. 
(8-16) Stellar mass [log(MQ)] with the associated uncertainty, estimated with the modeling parameters, [BC03,CHAB,CAL01] , and z-fit. 
(9-17) Stellar mass [logfM©)] with the associated uncertainty, estimated with the modeling parameters, [M05,KROU,CAL01] , and z-fit. 
(10-18) Stellar mass [log(Mo)] with the associated uncertainty, estimated with the modeling parameters, [CB09,CHAB,CAL01], and z-fit. 
(11-19) Stellar mass [log(Mo)] with the associated uncertainty, estimated with the modeling parameters, [CB09,SALP,CAL01], and z-fit. 
(12-20) Stellar mass [log(MQ)] with the associated uncertainty, estimated with the modeling parameters, [P01,SALP,CF00], and z-fit. 
(This table is available in its entirety in a machine-readable in the online version. A portion is shown here for guidance.) 



TABLE 9 

The IRAC-3.6+4.5 (mm sample: IR-Luminosities and Star Formation 
Rate estimates 



Object 


a 


<5 


f(24) 


f(70) 


z-fit 


SFR(0.28) 


SFR(R09) 


L(TIR,24) 


SFR(TIR.24) 


SFR(TIR,24,CE01) 


SFR(B08,24) 


SFR(AH06.24) 








err-f 


err-f 








L(TIR,best) 


SFR(TIR,best) 


SFR(TIR,CE01) 


SFR(B08,best) 


SFR(AH06,best) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 








(14) 


(15) 








(16) 


(17) 


(18) 


(19) 


(20) 


irac003270.1 


215.43892696 


53.08455063 


55 




1.09 


6.0 


21.8 


10.93 


14.6 


15.1 


17.5 


9.4 








13 










10.93 


14.8 


15.1 


18.4 


11.0 


irac003278 


215.42614011 


53.09447161 






0.49 


2.0 


1.7 


-10.34 


-3.8 


-4.7 


-8.1 


-2.2 


irac003291.1 


215.44043562 


53.08128671 


104 




0.86 


2.5 


17.4 


10.92 


14.4 


14.6 


17.5 


9.3 








12 










10.94 


14.9 


14.6 


18.5 


10.7 


irac003310 


215.42129738 


53.09430607 


61 




0.08 


0.1 


0.0 


8.61 


0.1 


0.1 


0.3 


0.1 








10 










8.51 


0.1 


0.0 


0.2 


0.1 


irac003313 


215.43553774 


53.08200958 






1.06 


3.6 


19.9 


-10.90 


-13.8 


-14.5 


-17.3 


-8.6 



Note. — (1) Object unique identifier in the catalog. 
(2,3) Right Ascension and Declination (J2000) in degrees. 

(4.14) Observed flux and uncertainty in MIPS 24 fj,m [/iJy]. 

(5.15) Observed flux and uncertainty in MIPS 70 fj,m [fijy]. 

(6) Rcdshift used in the fitting procedure. This rcdshift refer to zphot-best unless the spectroscopic rcdshift is available; in that case the r edshift is forced to the spectroscopic value. 

(7) UV based SFR [Moyr -1 ] estimated from the monochromatic luminosity at 2S00A rest-frame using the calibration of iKennicuttl (|T998).Wc also refer to this value as SFRt/v,o6s 

(8) IR based SFRs [MQyr - ] estimated from the observed flux in MIPS 24/im and the redshift using the formula o f lRieke et all i2009fl . 

("9,16) Total I R luminosity [log(L(? ) ' ) 1 obtained integrating (from 8-1000/xm) the average of the best fitting templates. (9) is computed from the fit of MIPS 24 fim data t o the models of | Charv fe Elbad 
d200ll CEOll. fDale & Heloul d2002l . DH02); (10) is computed from the fit of IRAC -8.0 and M I PS 24 and 70 /im (best effort) data to the models of CE01, DH02 and lRieke et all d2009l . R091. 

(10.17) IR based SFRs [MQyr -1 ] estimated from L(TIR) using the calibration of Kcnnicutt (1998). (10) and (17) are computed using the same combination of data and models as (9) and (16), 
respectively. 

(11.18) Same as 10, but in this case the IR SED is fitted only with the models of CE01. (11) and (18) arc computed from the fit of these modclcs to MIPS-24 fim data only and IRAC-8.0 and 
MIPS 24 and 70 f-tm respectively. 

(12.19) IR based SFRs [M yr _1 )] estimated from the monochromatic luminosity at Sfim rest-frame using the relation of Bavouzct ct al. (2008). (12) and (19) are computed using the same 
combination of data and models as (9) and (16), respectively. 

(13.20) IR based SFRs [Moyr - 1 ] estimated from the monochromatic luminosity at 24/j,m rest-frame using the realtion of lAronso-Hcr rcro et al . (2006). (13) and (20) arc computed using the same 
combination of data and models as (9) and (16), respectively. 

The uncertainties in the values of L(TIR) and the SFRs can be as high as a factor of 2. The accuracy in (9,16) is limited to two decimal places and one decimal place in (10-13,17-20). Negative 
values in the columns (9-13) indicate that the sources are non-detected in MIPS 24/im. In these cases, the corresponding IR luminosity is estimated from a upper limit of f(24)— 60/xJy, and the 
estimates in (16-20) arc not computed. 

(This table is available in its entirety in a machine-readable in the online version. A portion is shown here for guidance.) 
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